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Some Harmonic Properties of an Oscillating 
Fabry-Perot Interferometer 


M. Gadsden* and H. M. Williams* 


(March 9, 1966) 


The transmission of a spectral line through an air spaced Fabry-Perot interferometer is considered. 
The harmonic components of the transmission when the spacing of the interferometer is oscillated 
are discussed and it is found that the components can be used in a servocontrol system to offset 


drifts in the optical thickness of the interferometer. 


observations of the night airglow is discussed. 


A possible use of such a controlled system in 


Key Words: Fabry-Perot interferometer, airglow photometry. 


1. Introduction 


In two recent papers, Mielenz et al. [1, 2], have dis- 
cussed Fabry-Perot interferometers which use plates 
whose separation can be controlled electrically. The 
spacers in these interferometers are made from a 
piezoelectric ceramic material, lead zirconate-titanate, 
and the plate separation is thus directly, and easily, 
controllable by the application of an electric potential 
across two faces of the material. Mielenz et al., (loc. 
cit.) present full details of the use of this material and 
discuss the advantages and some applications of this 
piezoelectric scanning technique; it is the purpose of 
this paper to point out that the harmonic output of 
such a system has some advantages in the servo- 
controlling of the spacing of an interferometer. 

The use of oscillating piezoelectric spacers to con- 
trol the spacing and parallelism of a Fabry-Perot inter- 
ferometer was first described by Ramsay [3]. In 
Ramsay’s system, two passes of the interferometer 
by white light are used, and equality and constancy of 
the path length in the interferometer along the two 
light beams is maintained by seeking the maximum 
transmission for white light. The servocontrol to 
seek this maximum involves a minute superposed 
oscillation to give an error signal. The servocontrol to 
be described in what follows is similar to Ramsay’s 
system except that the control for constancy of inter- 
ferometer spacing is achieved from the spectral emis- 
sion being observed with the interferometer, and does 
not demand subsidiary white light illumination. 


*Institute of Telecommunication Science and Aeronomy, Environmental Science Serv- 


ices Administration, Boulder Colo. 


Formerly the Central Radio Propagation Laboratory 
of the National Bureau of Standards. 


2. Harmonic Analysis of the Output 


As is well known [4], the transmission at normal 
incidence of a perfect Fabry-Perot interferometer is 
given by the Airy relation 

: = ] 
Im 1+F sin? (27D) (1) 





where D= optical thickness, in wavelengths, of the 
interferometer i.e., half the order, 
/=transmitted intensity of a spectral line, 
whose width is much less than the instru- 
mental width: 


: ; n , 
I =transmitted intensity when D=Z=, n being 
_ any integer, 2 i 
F =a constant, dependent on the reflectivity of 
the coatings on the interferometer plates. 


F is related to the finesse, ¥, of the interferometer by 
the relation [4| 
~ (sry 
r=(*). 
7 


If a sinusoidal potential be applied to the piezoelec- 
tric spacers of an interferometer, D in expression (1) 
will become (dy + d cos 27 ft), where fis the frequency 
of the alternating potential and d is half the total 
excursion in length of the spacers. In these circum- 
stances, //I/m becomes a periodic function of time, f¢, 
with a period equal to 1/f. The function, //Jm, will 
contain a series of harmonics of f; we shall direct our 
attention to the amplitudes of the zeroth, first and sec- 
ond harmonics in what follows, with the convention 
that the nth harmonic has the frequency nf. 


and 
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The harmonic amplitudes and phases contained in 
I/Im will vary as 2d is an integer or not. This is 
readily seen from figure 1, in which is sketched, from 
top to bottom, 

(i) The variation of //Jm as a function of plate 
separation, do; 

(ii) two spacer oscillations, (A) and (B), with d= 1/2 

and d) at the positions indicated; 

(iii) the outputs resulting from the applied spacer 

oscillations, and 

(iv) two more outputs resulting from spacer oscil- 

lations with d=1/2 and do at the marked posi- 
tions, (C) and (D). 

It is readily seen that the second harmonic predomi- 
nates in (A) and (B) and that the first harmonic pre- 
dominates in (C) and (D). Photographs of oscilloscope 
traces showing approximations to these conditions are 
shown in figure 2; note, in both figure 1 and figure 2, 
that there is a change in phase of the output between 


(C) and (D). 


TRANSMISSION 





| 
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INTENSITY 








(C) 


FicuRE 1. Top: the variation in transmission of monochromatic 
radiation by a Fabry-Perot interferometer, as the separation of 
the plates is increased. 


Middle: the time variation of the transmitted intensity when the 
plate separation is oscillated over the ranges indicated by (A) 
and (B). 


Bottom: the time variations resulting from spacer oscillations 
equal in amplitude to (A) and (B), but centered at the points (C) 
and (D) indicated in the top sketch. 


FIGURE 2. Oscillograms of an oscillating interferometer corre- 
sponding to three settings of the mean spacing. 


Numerical harmonic analysis of //Im has given the 
results shown in figure 3. The amplitudes of the 
harmonics depend on whether 2d, is an integer, i.e., 
if we write 


dy=} (n + €), 


where € <1 and n is integral, all possible variations 


of do are covered by varying €. In figure 3, therefore, 
the abscissae are given in terms of €/2, with the limits 
0, + 0.5. 

The upper plot in figure 3 is for an interferometer 
with a finesse of 10, and d equal to 0.25. With this 
amount of driving potential, the point of interest is that 
the second harmonic is a maximum at the value of «€ 
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FIGURE 3. Harmonic amplitudes of transmitted monochromatic 
radiation as a function of the departure of the mean spacing from 
an integral of half-wave lengths. 


where the first harmonic is zero (and where, also, it 
changes phase). This, because of the symmetry of 
I/Im, holds for all values of the finesse so long as 
d=0.25. When d is either increased or decreased, 
the null point of the first harmonic remains at €/2=0.25, 
but the second harmonic shows two maxima at progres- 
sively increasing distances from €/2=0.25. (When d 
is increased to 0.5, then the two maxima of the second 
harmonic recombine at €/2=0, 0.5, where the first 
harmonic has another null point.) The lower plot in 
figure 3 shows the effect of reducing d to 0.15. 
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This property of the second harmonic going to a 
maximum where the first harmonic is zero and chang- 
ing phase, is of immediate application in using a 
servosystem to control the mean optical thickness, do, 
of the interferometer. dy is related to the metric 
thickness, d*, of the air layer between the interferom- 
eter plates by 


pd* 


do= hn ’ 


(3) 


in which vp is the refractive index of the air in the inter- 
ferometer and Ao is the wavelength in vacuo of the 
observed emission line. yp is, of course, dependent 
on the temperature and pressure of the air and, in 
the course of time, can be expected to show small 
changes; more important, however, is the temperature 
dependence of d* through the thermal expansion of 
the spacer material. It is easy, however, to compen- 
sate for these changes in d* and yp by arranging for 
a bias voltage to be applied to the spacer material 
from a servosystem which seeks the null of the first 
harmonic. With a tuned system responding to the 
second harmonic, one then has a controlled interfero- 
metric photometer for monitoring the radiance in a 
spectral line. 

In this connection, it is intuitively obvious that there 
will be a particular value for the finesse that gives 
the most output in the second harmonic. If the finesse 
is very high, the proportion of the cycle during which 
the spectral line is passed by the interferometer will 
be small, relatively more of the available signal 
As the finesse 


appearing in the higher harmonics. 
is decreased, the output signal will tend to become 
more and more sinusoidal, although with a finesse 
equal to zero—i.e., with no interferometer —there will 


be no second harmonic. Somewhere, therefore, there 
will be a maximum in the curve of the second harmonic 
amplitude versus finesse; figure 4 shows this maximum 
to be close to a finesse of six. 

In the above analysis, it has been assumed that the 
interferometer is ‘perfect, with J/Im given by (1). 
In practice, the performance of an interferometer 
falls short of perfection due to surface imperfections 
in the plates, lack of perfect adjustment in parallelism, 
and the use of a finite scanning aperture in the plane 
of the fringes. Measurements we have made using a 
prototype system, with a finesse equal to approximately 
five, have shown that the amplitude and phase rela- 
tionships between the harmonics, derived above, hold 
for a practical system. 

It-has also been found, during these measurements, 
that the ceramic material shows a ‘damped’ response 
to the application of a potential, such that if a poten- 
tial be switched on, the majority of the strain appears 
very rapidly but the final strain appears only after 
the lapse of a few seconds. This implies that the volt- 
age excursion required to produce a given excursion 
in length of the spacers will be a function of the fre- 
quency of an applied alternating potential, and will 
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FIGURE 4. Change of the maximum of the second harmonic with 
change in the finesse of the interferometer. 


be somewhat greater than that indicated from quasi- 

. static measurements. Since the first harmonic null 
point is independent of the amplitude of the alternating 
driving potential, a practical procedure of setting the 
bias and driving potentials is to set first the bias poten- 
tial to give zero first harmonic and then to vary the 
amplitude of the superposed alternating potential 
to maximize the second harmonic. We have found, 
using our prototype system, that this procedure is 
quick, and reliable. 


3. Airglow Photometry 


The problem basic to photometric determination 
of the night airglow is that of observing an atomic 
emission line spectrally embedded in a continuum 
and emission bands. Two basic techniques are in 
use to eliminate, or take account of, the radiation spec- 
trally surrounding an emission line. The most fre- 
quently used is that of ‘two-filter’ photometry, in which 
interference filters are used either to measure the 
radiance of the continuum in a spectral region free of 
atomic lines so that the contamination of the atomic 
line measurement may be estimated (Roach and 
Barbier [5]) or, in a method due to Barbier [6], to make 
two measurements centered on the emission line with 
differing proportions of background contamination. 
The first of these methods suffers from the uncertainty 
in knowledge of the spectral distribution of the back- 
ground radiation (Truttse [7, 8]); both ‘two-filter’ 
methods use pairs of observations sequential in time, 


and are thus inherently unsuited for use in observa- 
tions involving dynamic scanning of the sky. In this 
case, one uses the second basic technique of extracting 
an atomic emission line from the continuum, which is 
to use either a birefringent filter (Koomen et al. [9]), 
or a rocking interference filter (Filosofo et al., [10)). 
The calibration of a birefringent-filter photometer 
presents some difficulties although these are not in- 
superable (Smith and’ Alexander [11]; Dandekar et al., 
[12}). 

Let us consider the basic problem. Whena photom- 
eter is pointed at the night sky, the response, R, 
consists of two parts 


R=R,+ Re, (4) 


where R,, the response to the background continuum, 
is given by 


Ro=W-h-S (5) 


(W = equivalent width of the filter, /, = specific spectral 
radiance of the continuum and S= spectral response 
of the photocell). Similarly, R., the response to the 
emission line, is given by 


Re=Te-le-S> (6) 


(T.=transmission of the filter at the wavelength of 
the emission line and /,=radiance of the line.) 

If we now take an oscillating interferometer as the 
filter in the photometer, then 7,, given by (6) above, 
becomes of the form //I, discussed above; particularly, 
the first harmonic in R, can be used to compensate 
for thermal changes in the spacer length and the 
second harmonic used to determine /,._ Under these 
conditions, the contamination of the measurement by 
R,, eq (5), can become extremely small: first, by making 
W very small, as is easily done with an air-spaced 
interferometer, and, second, because the second har- 
monic in R, will depend only on the derivatives of 
(I,:S), taken at the emission wavelength, of even 
degree. These derivatives will be appreciable only 
if the emission line lies in a fraunhofer line in the 
continuum, or is close to emission band structure. 

This result is seen by expanding J, in Taylor’s 
series, ignoring the variations of S with A: 


Tio + AX cos wt)= Ip(Ao) + AA cos wt (S)a=% 


dX 


— dl, 
+ Ad? cos? wt (SE)a-r0+ “+ 


where (Ag+ AA cos wt) expresses the oscillating wave- 
length of the transmission peak of the interferometer 
pass-band. (It is to be noted that the coefficient of 


d?1), dl], 


cos 2 wt contains onl =. ier 
y dd?’ dd4 


-- etc.) 
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It should be noted that the presence of appreciable 
odd derivatives in (J,-S) will give a contribution to 
the first harmonic signal, causing the servosystem to 
seek a null point away from the peak of the second 
harmonic. In the presence of a serious background 
continuum, therefore, it would be best to make 
W as small as convenient. 


4. Conclusions 


A system such as the one outlined above offers 
the possibility of continual control of the spacing, 
while observations of a feeble spectral line are in pro- 
gress. The use of three small spectral line sources, 
sampling the interferometer at three separated posi- 
tions, is an obvious extension of this monitoring to the 
control of parallelism. In this case, we have a system 
closely resembling that of Ramsay [3], but one which 
does not demand an optical system, with close toler- 
ances, to return the sampling light for a second pass 
through the interferometer. 
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Precise Reflection Coefficient Measurements With an 
Untuned Reflectometer 


W. E. Little* and D. A. Ellerbruch* 


Institute for Basic Standards, National Bureau of Standards, Boulder, Colo. 


(March 9, 1966) 


The precision tuned reflectometer technique of reflection coefficient measurement is, at present, 


the most accurate technique available. 


Utilization of the technique requires the use of tuners and 


sliding terminations to eliminate or reduce the inherent errors of a directional coupler. 
This paper describes a reflection coefficient measurement technique that approaches the accuracy 


of the tuned reflectometer. 


Instead of tuners, the technique uses an ideal one-quarter wavelength 


waveguide section to eliminate the measurement error introduced by finite directional coupler direc- 
tivity and to reduce to second order the error introduced by an equivalent generator mismatch. 

A comparison of the quarter-wave technique with the tuned reflectometer technique through the 
measurement of the reflection coefficient of three 9/16 inch coaxial waveguide terminations is also 


included. 


7 . bor 
Key Words: Microwave, impedance, measurements, reflectometer. 


1. Introduction 


The reflectometer technique for microwave reflec- 
tion coefficient measurements is now widely used 
[1]... The measurement accuracy of a basic reflectom- 
eter (a simple directional coupler coupled to the re- 
flected wave) depends on the directivity of the 
directional coupler and the equivalent generator mis- 
match. The best known method of improving the 
measurement accuracy of the basic reflectometer 
requires the use of tuners, tuning procedures, and 
auxiliary equipment [2, 3]. The result of this tuning 
is an increase in the effective directivity of the coupler 
and a matched equivalent generator. At present, 
the most accurate measurements are achieved with 
this tuned reflectometer. 

The purpose of this paper is to present another 
method of improving the measurement accuracy of 
the basic reflectometer technique. This method 
utilizes a quarter-wavelength . waveguide section. 
Proper use of the quarter-wave section during the 
measurement procedure will, in theory, eliminate 
the error caused by a finite directivity and will reduce 
the error caused by a mismatched equivalent generator 
to second order. 

No detailed error analysis is included in this paper; 
however, a laboratory measurement comparison with 
the tuned reflectometer indicates that the total meas- 
urement error is small. Such an error analysis may 
be given at a later date, depending on results of further 
laboratory experiments with the method. 


* Radio Standards Engineering Division. National Bureau of Standards, Boulder, Colo. 


' Figures in brackets indicate the literature references at the end of this paper. 


2. Theory 


The operation of a reflectometer is discussed in 
detail elsewhere [3]. Therefore, only a brief review 
will be presented here as a starting point for the dis- 
cussion of the quarter wavelength technique. The 
side-arm output, bs, of a directional coupler connected 
to sample the reflected wave (see fig. 1) can be ex- 
pressed as 


— 
gtlr 


bs = kbe 1—faly’ 


(1) 


where b; is the component of the incident wave fur- 
nished by the generator, I'2; and k are both functions of 
the scattering coefficients of the directional coupler 
and reflection coefficients of the generator and detec- 
tor, and Iz is the reflection coefficient of the terminat- 
ingelement. The symbol I2; represents the equivalent 
source reflection coefficient presented to the reflector 
under study at the reference plane. The symbol K 
is approximately equal to the directivity ratio of the 
directional coupler. 

The desired performance of the directional coupler 
is realized when I,;=0 and 1/K=0. Equation (1) 
then reduces to 

bs = kb<Ir. (2) 
Hence, inspection of (1) shows that the terms involving 
1/K and I; introduce first order errors in 7, when Ir 
is determined according to (2) from a measurement of 
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the side-arm output, 63, of an untuned (1/K and T2; 4 0) 
directional coupler. Reference 3 discusses the mini- 
mization of these errors by proper use of tuners. 


It is possible to eliminate the term 1/K and simul- 
taneously reduce the error arising from a finite ['2; to 
second order by a measurément procedure utilizing an 
ideal waveguide section that is one-quarter wavelength 
long. The theoretical basis for the method is as 
follows: 

If the unknown is connected directly to the reference 
plane of the untuned reflectometer shown in figure 1, 
the side-arm output, bj, is 


Lyn, 


bu = kbc TT, 


(3) 


where I’, is the reflection coefficient of the unknown. 

Interposing a precision quarter-wavelength wave- 
guide section between the reflectometer and I, adds 
a 180° phase shift to the argument of the impedance 
appended to the reference plane, causing the reflec- 
tion coefficient at this plane to be —I,. With the 
quarter-wavelength section inserted between the refer- 
ence plane and the unknown, as shown in figure 2, the 
side-arm output, bo, is 








=! 
| 
porns 


REFERENCE 
PLANE 


FIGURE 1. A basic reflectometer. 








REFERENCE 
PLANE 
FIGURE 2. 


A basic reflectometer with a 


connected. 


quarter-wave section 


The difference between b,, and bz, is then 


- 
I+ Psi 


bi, — bo, = 2kb. ly 1-7," 
24 pl)” 


(5) 


If the unknown is now replaced with a quarter- 
wavelength standard short circuit? whose reflection 
coefficient is unity and the above measurement 
procedure repeated, formula (5), upon substitution 
of [;=1 for T,, reduces to 


+2. 


bis — bos = 2hby; (6) 


1 —(P2;)?’ 


where b,; is the sidearm output when the standard 
short circuit is connected directly to the reference 
plane, and b.; is the sidearm output when the ideal 
quarter-wavelength section is placed between the 
reference plane and the short circuit. 

Dividing (6) by (5) results in 


bis — bes = 5 ] — (TaiJ0°,,)? . 
kicwtes ty tke 





(7) 


Note that (7) does not contain the 1/K factor and that 
I'2; occurs only to the second power. |I2;| is usually 
less than 0.05 in a practical situation; hence, | I>; |? 
<1. Therefore, the limitations of the untuned direc- 
tional coupler (i.e., finite directivity and a mismatched 
equivalent generator) have been essentially circum- 
vented, and the correct value of I’, can be determined 
without tuning the coupler. 

To evaluate [, by means of (7), it appears that four 
individual (complex) quantities have to be measured. 
If one is interested only in a magnitude measurement, 
which is usually the case, a relatively simple measure- 
ment procedure can be used to determine the ratio 
| bis — bes |/| bi,—be,| directly, in decibels. 


3. Measurement Procedures for Determining 
the Modulus of I’, 


A measurement system has been constructed and 
used to substantiate the quarter-wave method experi- 
mentally. A block diagram of the system is shown 
in figure 3. The measurement procedure is as 
follows: 

1. One end of the precision quarter-wavelength 
section of waveguide is connected to the reference 
plane and terminated with a standard short circuit. 
The signal from the side arm under these conditions 
is bo, of formula (7). 


2A quarter-wavelength standard short circuit is a quarter wavelength of precision line 
terminated by a shorting metal wall. If a simple shorting plate is used, [,=— 1 is sub- 
stituted for [’, in (5) and one side of (7) would be multiplied by a — 1. 
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FIGURE 3. Block diagram of quarter-wavelength measurement 


system, 


2. The variable phase shifter and attenuator, 
designated P, and A, in figure 3, are adjusted to 
produce a null at the 30 MHz IF receiver. This 
adjustment sets the signal from the auxiliary channel, 
bia, equal to — bos. 


3. The ideal quarter-wavelength section is removed 
and the standard short circuit is connected directly 
to the reference plane. Under these conditions the 
signal from the side arm is bis of (7). Because the 
signal from the auxiliary channel of figure 3 was ad- 
justed in the preceding step to produce the signal 
—b»;, the receiver, which responds to the sum of the 
outputs of the side-arm and auxiliary channel, registers 


| bu bos|. 


4. The standard 30 MHz attenuator in the IF re- 
ceiver is adjusted to produce an arbitrary reference 
level at the indicator. 


5. The short circuit is removed from the reference 
plane, and the quarter-wavelength section is connected 
and terminated in the unknown reflection coefficient. 
Under these conditions the signal out of the side- 
arm is equal to bo, of (7). 


6. The phase shifter and attenuator (P; and A;) are 
now readjusted to again produce a null at the 30 MHz 
receiver. This adjustment provides that b4.=— boy, 
where 6,» is the signal from the auxiliary channel. 


7. The quarter-wavelength section is removed, and 
the unknown is connected directly to the reference 
plane. The signal into the 30 MHz IF receiver is 
again the sum of the signals from the side-arm and 
auxiliary channel, i.e., |bi,—62,|._ A record is made 
of the change (in decibels) in the setting of the stand- 
ard 30 MHz attenuator in the receiver required to 
return the indicator to the reference level chosen in 
step four. This difference (in decibels) between the 
final attenuator setting and the setting at step four is a 
measure of the ratio |b1;— bes|/|bi,—b2,|. Thus, 


dB(measured) = 20 log ete 
Iw” D2 


= mT 11 —(P2T,,)?| 
ler] T—Bi 


Assuming that 


this equation reduces to 
approximately 


dB(measured) = 20 log =s 


Pal’ 


whence the unknown reflection coefficient magnitude 
is given by 


__ dB(measured) 
20 ° 





[P.| = 10 


4. Measurement Results 


The reflection coefficient magnitude of three 9/16- 
inch coaxial terminations were measured at 4 GHz 
using the quarter-wave technique. These measure- 
ments were compared with those obtained with a 
tuned reflectometer and those obtained with an un- 
tuned reflectometer that did not utilize a quarter- 
wavelength section. Table 1 summarizes the results. 

The first column lists the values obtained for the 
three terminations when measured with a _ tuned 
reflectometer. The tolerances associated with these 
values include the measurement error as well as the 
non-repeatability obtained by rotating the termination. 
In the second column are listed the VSWR’s corre- 


sponding to the reflection coefficients in the first 
column. 


TABLE 1. Tabulated measurement results 





VALUES OBTAINED WITH 
IMPERFECT: COUPLER 
AND !142 SECTION 


IT | 


TUNED REFLECTOMETER 
MEASUREMENT VALUES 


IT! 


VALUES OBTAINED WITH 
IMPERFECT COUPLER 


ITI 





VSWR VSWR VSWR 





0.0217 
+0.0006 


1,0443 
+0.0012 


0.042 1,088 0.0224 1.0458 


0.0902 
+0,0015 


1.1983 
+0.003 


0,108 1.242 0.0918 1,2022 


0.3334 
+0,0028 


2.000 
+0.006 


0.351 2.082 0.3353 2.0089 


























In the third and fourth columns are listed the re- 
flection coefficients and VSWR’s obtained with an 
untuned directional coupler but no quarter-wavelength 
section. These values were obtained as follows: 
The standard short circuit was connected to the 
reference plane, and an indicator reference point 
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was obtained. The short circuit was then replaced 
by the unknown impedance, and the reflection co- 
efficient was computed from the measured change in 
the signal level at the indicator. These values are 
subject to all the errors due to imperfections in the 
directional coupler. The directional coupler employed 
in the three systems with which the data of table 1 
were obtained (including that of fig. 3), had a direc- 
tivity of approximately 30 dB (this corresponds to a 
1/|K| value of approximately 0.032) and a |Ts;| of less 
than 0.02. 

In the fifth and sixth columns of table 1 are listed 
the reflection coefficient magnitude and VSWR values 
obtained by using the quarter-wavelength measure- 
ment technique described in this report. 

The values obtained using the quarter-wave tech- 
nique are in near agreement with the values obtained 
using the tuned reflectometer. The values obtained 
with the untuned reflectometer when the quarter- 
wave section was not used are considerably in error. 


5. Conclusions 


Precise reflection coefficient measurements are 
possible with an untuned reflectometer if the quarter- 
wave measurement technique is used. The theory 
was verified experimentally by measuring the reflector 
coefficient magnitudes of several terminations with a 
tuned reflectometer and then remeasuring them 


using the quarter-wave technique with an untuned 


reflectometer. These same terminations were also 
measured directly with an untuned reflectometer with- 
out application of the quarter-wavelength technique, 
to demonstrate the improvement in accuracy provided 
by the quarter-wave method. 

Relative phase data were neglected during this 
evaluation. No serious difficulties are anticipated 
in modifying the measurement system to obtain phase 
data so that the quarter-wave technique will then 
provide the complex reflection coefficient. 

Use of the quarter-wave technique instead of a tuning 
procedure eliminates the need for variable-phase 
high- and low-reflection loads and shortens the required 
length of precision line. The most important potential 
application of the quarter-wave technique is in the 
measurement of impedance at frequencies below ap- 
proximately 1 GHz. With decreasing frequency the 
tuned reflectometer becomes increasingly cumber- 
some. 
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A 9-dial resistance ratio instrument capable of precise high ratio comparison of voltage or re- 
sistance (from 1/1 to 107/1) is described. Also described is the modification of a 6-dial universal 


ratio set to permit its additional use as a versatile ratio instrument. 


Paramount to the accuracy of 


these high ratio instruments is the carefully adjusted “common point” junction which is briefly 


discussed. 


Key Words: DRRS (direct reading ratio set); junction; precision measurements; ratio, resistor, 
four-terminal; resistance decade; transresistance; URS (universal ratio set); VERI 


(versatile ratio instrument). 


1. Introduction 


Because of a need to make accurate high ratio d-c 
resistance comparisons and an inability to do so with 
existing ratio instruments such as the DRRS (direct 
reading ratio set) and the URS (universal ratio set), 
[1, 2, 3]! a versatile ratio instrument (abbreviated VERI) 
was developed. With the VERI an arbitrary voltage 
or resistance ratio from 1/1 to 107/1 (or 1/107) may be 
measured with a dial resolution of 1 part in 108 or 
better. (Permissible generator voltage and detector 
sensitivity may further limit the actual resolution at- 
tainable.) This instrument may be used to calibrate 
ratio devices — potentiometers, voltage dividers, volt 
boxes, the ratio arms of Wheatstone bridges, etc. 
Also, the individual fixed resistors of the VERI may 
be used as 4-terminal resistance standards, and the 
incremental resistance section may be used as a 9- 
dial, 2- or 4-terminal resistance decade. 

A 6-dial URS may be modified as described in sec- 
tion 3 to permit its additional use as a 6-dial VERI 
having a ratio range from 1/1 to 2 X 104/1. Its normal 
function as a URS is not impaired by the modification. 

The accuracy of the VERI depends directly upon 
the design and adjustment of the “common point” 
junction [4]. The transresistances? of the junctions 
used in the described instruments are less than 10-8 2. 
This causes a worst-case uncertainty in ratio of less 
than 1 part in 10.’ 


*Radio Standards Engineering Division, National Bureau of Standards, Boulder, Colo. 


1 Figures in brackets indicate the literature references at the end of this paper. 
2 Transresistance: The ratio of the voltage induced at one pair of junction terminals to 


the current through another terminal pair, i.e., Rue=73, see figure 1. 
13 


JUNCTION TO BE MEASURED 
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FIGURE 1. Basic setup for measuring junction transresistance. 


2. Versatile Ratio Instrument 


Figure 2 is a schematic diagram of the VERI. It 
consists basically of a precision resistance decade 
section and several precision fixed resistors of decimal 
values connected at a common point by a “zero” 
transresistance junction. The resistance decade sec- 
tion shown consists of 9 decades which may be varied 
in 0.001 © steps to a total resistance of 1,111,111.111 
Q. The seven fixed resistors have decimal values 
from 0.1 to 100 kM. The “outside” ends of each 
of the resistors have dual terminals, so they may 
be connected as 4-terminal resistors to permit the 
elimination of the terminal and external lead resist- 
ances. The “common-point” junction shown is a 
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FIGURE 2. Schematic diagram 
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composite junction consisting of three tetrajunctions,* 
A,B,andC. This is necessary because a junction may 
have a maximum of only four electrically symmetrical 
terminals [4]. Junction A is the base junction to 
which all resistance and ratio measurements are re- 
ferred. It is adjusted to have less than 10-* 2 
transresistance. Because of their relatively large 
resistance, the 10 kN and 100 kO fixed resistors, and 
the 0-1 MQ resistance decade section may be con- 
nected to a common terminal of the base juction 
without introducing significant error (assuming ra, 
is less than 0.001 Q). Terminals a, b, and ¢ are 
needed for the calibration and use of the VERI since 
the potential terminal used for a ratio measurement 
should not be common to either of the resistors being 
used in the VERI. The switches used in the resist- 
ance decades have low resistance variation (less than 
0.001 2 per switch). 

This instrument when fully calibrated may be used 
in a variety of ways, including: 

(1) 4-terminal fixed resistance standards, 0.1, 1, 
10, 100, 1k, 10k, and 100 kQ: 

(2) a 2- or 4-terminal, 9-dial decade resistor, 0.000 
to 1, 111, 111.111 0; 

(3) a 0.01 ppm resolution direct reading ratio set for 
precise 1 to 1 comparisons; 

(4) fixed ratio standards with nominal ratios to 10®/1; 
and 

(5) an incremental ratio standard having continuous 
ratio adjustability from 1/1 to 107/1 with 9-dial 
resolution. 

Calibration of the instrument may be accomplished 
by measuring one of the fixed resistors, and then 
making 10/1 [5] and/or 100/1 step-up or step-down 
(ratio) calibrations of the remaining fixed resistors 
using the technique reported by Hamon [6| and further 
analyzed by Page [7] or, each of the fixed resistors 
may be measured and the ratios between them cal- 


*A “tetrajunction” is an electrically symmetrical.4terminal junction [4] having essentially 
zero transresistance. 


culated. The latter has the disadvantage that the 


uncertainty of the calculated ratios will include the 
uncertainty in the values of the individual resistors. 
The decade resistors are calibrated as a 4-terminal 
resistance box with direct checks of the xl and x10 
dial settings against the fixed resistors. 


3. URS Modified to a VERI 


A universal ratio set may be modified to permit its 
additional use as a versatile ratio instrument. A minor 
modification of a URS for use as a DRRS has been 
described by Ramaley and Shafer [8]. Figure 3 is a 
photograph of a URS which has been altered as shown 
by the schematic diagram of figure 4. As modified, it 
may be used as: 

(1) a conventional 6-dial URS: 

(2) 4-terminal fixed resistance standards, 0.1, 1, 10, 
100, and 1000 2: 

(3) a 2- or 4-terminal 6-dial decade resistor, 0.001, 
to 2111.110 O; 

(4) a 1 ppm resolution direct reading ratio set for 1/1 
comparisons; 

”) fixed ratio standards with nominal ratios to 104/1; 
an 

(6) an incremental ratio standard having continuous 
ratio adjustability from 1/1 to 2x 104/1 with 6-dial res- 
olution. 

The modifications to the URS consist of the following: 

(1) Replacing the x100 2 switch (top decade) with 
a switch having very low contact resistance varia- 
tion (less than 0.001 )). 

(2) Adding 12 terminals to the front panel, making 
a total of 17 as seen in figure 3. 

(3) Mounting a tetrajunction (a /2-inch diam copper 
sphere having four terminals). 

(4) Mounting the five additional precision resistors, 
0.1, 1, 10, 100, and 1000 2. 


(5) Wiring the resistors and terminals as shown. in 
figure 4. 
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FIGURE 4. Schematic diagram of the NBS-modified 6-dial URS to 


permit its additional use as a VERI. 


The spherical tetrajunction, the five added pre- 
cision resistors, and some of the additional wiring 
can be seen in the photograph of figure 5. 


4. A Typical Application of the VERI 


Figure 6 shows a modified 6-dial URS used as a 
VERI in a setup for calibrating the ratio arms of a 
Wheatstone bridge. Units L,; and Lz are lead com- 
pensating networks which are needed when resist- 
ances of 1000 or less are involved in the ratios. S, 
and S, are low resistance shorting wires which are 
used during the adjustment of the lead compensators, 


17] 


a t.35 34% 


A commercial 6-dial URS which has been modified to permit its additional use as a VERI. 


L, and L,. Balance is achieved and the setting of the 
R, dials of the URS is noted when the removal and 
replacement of the shorting links, S; and S2, have 
no effect on the galvanometer null. The ratio of the 


Wheatstone bridge ratio arm for the circuit values 
shown is 


_Ralt+a)_ Ra 


M=o1a+b)~ 0.1 


(1+a—b), 


where R, is the “R,’ reading of the URS dials, and 
a and b are the proportional corrections to Ra and 
0.1 Q respectively. (The right-hand expression of 
the equation introduces less than | part in 10° uncer- 
tainty when a and 6 are each less than 0.0001.) 

Once the ratio arms of the Wheatstone bridge are 
calibrated, the standard arm, Rs, may be calibrated 
by measuring the fixed and decade resistances of 
the VERI on the “unknown” terminals of the Wheat- 
stone bridge. 


5. Conclusion 


The versatile ratio instrument described permits 
the precise comparison of ratio from 1/1 to 107/1 
with a dial resolution of 1 part in 108 or better. By 
comparison, a standard 6-dial universal ratio set 
cannot measure a ratio of 107/1, and will give only 
3-digit resolution on a ratio of 104/1. A VERI modi- 
fied URS, however, will give 7-digit resolution on a 
ratio of 104/1 and not less than 6-digit resolution on 
any ratio from 1/1 to 2 x 104/1. The effective resolu- 
tion attainable in high ratio measurements may be 
limited by the maximum permissible generator (bat- 
tery) voltage and detector (galvanometer) sensitivity. 
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FIGURE 5. Internal view of the commerical 6-dial URS modified by NBS to permit its additional use as a VERI. 
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FIGURE 6. Circuit used to calibrate the ratio arms of a Wheatstone 


bridge with the VERI (modified URS). (Paper 70C3-228) 
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The comparison of volt boxes within uncertainties of +0.01 percent or less has, in the past, been a 


long and difficult task at the NBS Boulder Laboratories. 


A console for this purpose recently put in 


use has reduced the time and labor to a great extent and, at the same time, has reduced the uncertain- 
ties associated with the measurement circuitry to less than 0.001 percent, thus allowing realizable 


calibration uncertainties of less than 0.005 percent. 
which the operator was considered to be the decision-making link in the measurement chain. 


The console was designed as part of a system in 


The 


circuitry, through the use of a contained computer and associated circuits, provides the operator with 


the data he requires in the form he chooses. 


The operator is thus freed from the need to consider 


minutiae, and thereby is able to concentrate on those factors requiring judgment. 


Key Words: Calibration, console, electrical, measurement, ratio, standards, voltage, volt box. 


1. Introduction 


A volt box is a carefully constructed voltage divider 
having a number of fixed ratios between input and 
output terminals. These ratios are commonly in 
multiples of the integers 1 through 5. A volt box is 
normally used to extend the effective range of a po- 
tentiometer so that voltages up to 750 or 1500 V may 
be measured. 

As a volt box is a ratio device, the absolute resistance 
of its elements is not of great importance and under 
ordinary circumstances need be known only in a 
nominal sense (i.e., 750 0/V, etc.). However, before 
a volt box may be used for accurate measurements, 
the actual ratio must be determined. There are two 
characteristics inherent in volt boxes that manifest 
themselves in a dependence of the ratio upon the 
applied voltage. One is the change in the resistance 
of the elements with heating. In the large “‘master” 
volt boxes fabricated in accordance with NBS RP1419 
[1],! the changes in ratio resulting from self-heating 
are usually quite small, being in the order of 5 to 20 
ppm (parts per million) at rated voltage. In the less 
complex, general purpose boxes, however, the differ- 
ence between the ratio at reduced voltage (under 30% 
of rated) and that at rated voltage may be as great as 
200 ppm. The other characteristic is the variation 
in the resistance of the insulating parts of the box. 
Even when using the new, high-resistance insulating 
materials, small currents across the surface of the 
insulator can result from surface contamination. 
When these currents are in parallel with resistance 
paths, the ratio is disturbed. These currents are not, 
in general, a linear function of applied voltage, and 


*Radio Standards Laboratory, National Bureau of Standards, Boulder, Colo. 
' Figures in brackets indicate the literature references on page 180. 


the current paths cannot be well-defined. To reduce 
the effects of these currents, the more complex, 
master boxes have guarded circuits wherein the 
voltage across the leakage path is reduced considerably 
with a corresponding reduction in leakage currents. 

When these characteristics are considered, it be- 
comes apparent that the ratio should be measured at 
the voltage for which the box will be used. This may 
be done most easily by comparing the box with another 
volt box having known characteristics. Two methods 
of comparison are described in RP1419. (The second 
method is discussed in appendix I of that paper.) The 
first method is the one used at NBS. 

The recently completed console used for this meas- 
urement at the NBS Boulder Laboratories incorpo- 
rates many new features that contribute to the accuracy 
and efficiency of the measurement. With this console, 
uncertainties associated with the measurement circui- 
try are less than 0.001 percent, which allows realizable 
calibration uncertainties of less than 0.005 percent. 
This paper describes the console and includes, as 
appendices, an evaluation of the uncertainties that may 
be expected in its use. 


2. Comparison Measurement 


As a comprehensive description of this measurement 
may be found in the literature, [2] only a brief treatment 
will be undertaken here. The circuit is shown in 
figure 1. In the course of the measurement, the values 
of v3, v2, v;, and vy are measured in microvolts. (vs is 
normally measured first because it is the voltage most 
sensitive to changes in ratio caused by heating and to 
wrong connections in the measuring circuit.) When 
these values and the supply voltage, V, are known, it 
can be shown [3] that the corrected ratio of the box 
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FIGURE 1. Basic volt box measurement circuit. 


undergoing comparison (within the limits of uncer- 
tainty in the measurement) is 





ry v4— U1 V2 — V3 6 
X Xu) 1+ (ot 7 +2) 10 |. 


where X,, is the nominal ratio of the unknown, 


fs is the correction to the ratio of the standard volt 
box in ppm, 

S» is the nominal ratio of the standard volt box, and 
is equal to Xp, 

V is the supply voltage, in volts, and 

V1, V2, V3, and vy are the measured voltages, in 
microvolts. 


(v;, v2, and vz are considered to be positive when of 


the same polarity as v3.) 


3. Measuring System Requirements 


In order to achieve a high degree of accuracy, the 
measuring system must, of course, be fabricated of 
stable, high-quality components and must be guarded. 
However, the use of good quality equipment alone is 
not a guarantee of highly accurate measurements. 
The success of the measurement is dependent to a 
very large extent upon the operator. Human opera- 
tional factors contributing to degradation of accuracy 
are primarily fatigue, both muscular and nervous, 
and inability to copy figures indefinitely without error. 
In the design of the console, these factors were care- 
fully considered, and as a result the following re- 
quirements were formulated: (1) The operator should 
be able to change ratios quickly and easily. This 
will not only reduce irritation but also will minimize 


the time during which no power is applied and thus 
reduce uncontrolled cooling of the volt boxes. It is 
important to minimize this cooling because volt boxes 
are calibrated under temperature equilibrium condi- 
tions. (2) Controls should be centralized as much 
as possible. (3) The digital form of indication is to 
be preferred over the analog form. (4) All readings 
and settings should be capable of being checked. 
(5) As much of the computing as possible should be 
done automatically. (In the calibration of a volt box 
eq (1) may have to be calculated as many as 32 times.) 
(6) Recopying of data should be reduced to a minimum. 

In addition to the requirements listed above, the 
fact that voltages as high as 1500 V dc may be present 
during calibration dictates the use of adequate safety 
measures. 


4. Features 


The console was designed from the human engi- 
neering viewpoint to as great an extent as was 'tech- 
nically and economically practical. The operator 
was treated as a decision maker, and the console was 
designed to automatically handle those functions not 
requiring decision and to supply the operator with the 
appropriate data in a clear and concise form (fig. 2). 
As a result of this viewpoint, and of the requirements 
listed above, the following features were included. 
(1) A plug-panel was designed to allow the operator 
to quickly change ratios on the standard volt box. 
(2) All controls are grouped in front of the operator. 
During the measurement of a given ratio, the operator 
does not need to leave his seat, and no control is more 
than 18 in. away. However, the operator does need 
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FiGURE 2. The volt box calibration console at the NBS Boulder 
Laboratories, now at NBS laboratories, Gaithersburg, Md. 
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FiGURE 3. A. Main control panel; B. Output 


to get up and move around the end of the console to 
change ratios. (3) All important control settings and 
displays are in digital form (fig. 3). (4) The basic 
measuring circuit includes provision for reducing the 
measured voltage to parts per million continuously 
and automatically. The value, in ppm, is recorded 
by the operator on an adding machine. (5) The tape 
from the adding machine, when complete, is cut into 
appropriate lengths and photographically reproduced 
on 842 X 11 inch sheets which become the permanent 
records. Thus the data need be recopied by hand 
only once, onto the summary sheet from which the 
Report of Calibration is typed. To comply with the 
safety requirements, the area in which the high voltage 
is exposed is protected by a two-stage photoelectric 
fence. If the illumination to either of two photocells 
is interrupted, the high voltage is disconnected. 
Also, a prominently mounted, high-voltage sign is 
illuminated whenever the high voltage circuit is ener- 
gized. To provide an indication of high voltage even 
if all other safety circuits were to fail, two neon bulbs 
are connected, in series with suitable dropping re- 
sistors, across the high voltage leads. One of these 
is mounted on the front panel, and the other is mounted 
at the top of the plug-panel in the high-voltage area. 
Initially a digital printer was connected to print the 
position of every control on the main control panel 
and the output from the measuring circuit. The tape 
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display panel. 


from the printer was to be used to check the validity 
of the recorded data. However, the action of the 
printer adversely affected the measuring circuit, and 
it proved to be impractical to reduce this circuit inter- 
action to a tolerable level. The use of the printer 
was therefore discontinued. The checking of re- 
corded data is accomplished by the making of multiple 
measurements. 


5. Measurements Procedure 


In the measurement of a given range of a volt box, 
the operator first connects the box as shown in figure 2. 
The main function switch (the large, rotary switch 
in the lower center of the main control panel, fig. 3A) 
is rotated to the “Enter S,”’ position. The operator 
enters the nominal value of the ratio into the console 
circuitry by adjusting the switches in the upper right 
of the main control panel. He then rotates the main 
function switch to the “Enter and Set V”’ position. 
(In so doing he moves the switch through the ‘Print 
S,”” position which originally caused the printer to 
record the values of the S, dials. This position no 
longer has any meaning.) After having turned the 
four ‘““V=’"’ dials to the desired value of applied voltage 
(the fourth dial controls the position of the decimal 
point), he adjusts the high-voltage supply until a null 
is obtained on the Set-V meter. By rotating the main 
function switch to the “Read v3” position, he makes 
those circuit connections necessary for the measure- 
ment of the difference voltage shown in figure 1. He 
records the values of V and S, on the calculator, and 
then proceeds to measure v3, then v2, v;, and vs. As 
the measurement of each difference voltage is com- 
pleted, he enters the value displayed on the output 
display panel (fig. 3B) into the calculator. After having 
completed the four measurements, he enters the cor- 
rection to the standard into the calculator. By de- 
pressing the totaling key, he completes the calculations. 


6. Computer Principle 


If eq (1) is rewritten as shown in eq (2), the correc- 
tion to the ratio of the box undergoing calibration is 
seen to be the sum of five terms, each of which is in 
parts per million, provided the conditions listed under 
eq (1) are met. 


= Us U1, Snv2 
X=Xy[1+(ue+$ yt 


=a) 10-8). (2) 


As only multiplication and division are needed to 
convert the measured voltage to parts per million, 
an electrical analog is easily constructed. Consider- 
ing the vs term of eq (2), 


S. 
Vx0ppm) = Va Ty (3) 
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Correspondingly, in a voltage divider circuit, 


R2 


Eo= ER: 


(4) 


where Rz is the portion of R; across which Eo appears. 
It remains only to determine the circuit that allows 
the relationship 


Eo _ R2_ Sn 
ER, V (9) 


to hold, considering that V and S, are independent 
parameters. 

If R» is allowed to vary in such a way as not to change 
the value of R,, the two values will be independent and 
eq (5) may be realized in fact. This is done as is 
explained in section 7.3. The result is a simple, 
completely passive computer, the accuracy of which 
is limited only by the accuracy of the contained 
resistors. 

In the actual design, no attempt was made to retain 
the decimal point in the voltage divider circuit, as 
that would require excessive voltages in the computer. 
Instead, a separate switching circuit is used to place 
the decimal point. 


7. Circuit 


7.1. Block Diagram 


Figure 4 shows a block diagram of the circuit. 
The output from the 0-1500 V supply is monitored 
by the set-volts potentiometer. This potentiometer 


is also used to enter the significant figures of the value 
of V (eq (2)) into the computing potentiometer such 
that R, of eq (5) is proportional to V. The signi- 
ficant figure of S, is similarly entered by adjusting the 
value of Re» in eq (5) appropriately. Decimal infor- 
mation is supplied to the decimal point computer as 
shown. (The value of V may be anywhere between 0 
and 1500 V, to three significant figures. That of S, 
may lie anywhere between 1.X10° and 9. xX 104, 
but in practice it does not exceed 1 X 104.) The four 
difference voltages (v;, ve, v3, and v4) to be measured 
are connected to the computing potentiometer by 
means of a relay switching #ircuit. The output from 
the computing potentiometer is connected to a digital 
volt-meter which has been modified to read only signi- 
ficant figures. The decimal information proceeds 
directly to the digital display from the decimal point 
computer. The printer monitored both the computing 
potentiometer (which included, for the purpose of 
printing, the settings on the set-volts and S, dials) 
and the digital display. The printer output included 
codes indicating the nature of the display (ppm, 
microvolts, etc.) as well as the value and sign. 


7.2. Set-Volts Potentiometer 


a. Significant Figures Section 


The set-volts potentiometer circuit is shown in 
figure 5. The significant figures section comprises 
a one-volt zener source supplying a Kelvin-Varley 
voltage divider. Each decade of the divider is 1000 2 
per step with 2500 (2 resistors used for compensation. 
This type of design allows all decades to be of a high 
(and easily procured) resistance which reduces the 
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FIGURE 4. Block diagram of console circuit. 
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effect of contact resistance. Each decade is physi- 
cally ganged to, but guarded and shielded from, a 
corresponding decade in the computing potentiometer. 
Thus, as the value of V is set on the three decades in 
the set-volts potentiometer, the value of R, in eq (5) 
is simultaneously adjusted to be proportional to V. 
It should be noted that the significant figures of the 
lowest value that can be set on the set-volts poten- 
tiometer dials are 100. This is necessitated by the 
need, in the computer, for Rigmin = 10 kQ > Ramax) 
=9kQ. A more detailed explanation is given in sec- 
tion 7.3. 

It will be noted that the galvanometer, G, (which is 
of the chopper-amplifier type), is protected by three 
fuses (fi, f2, and fs) in series. The possibility, how- 
ever remote, of the operator inadvertently applying 
1500 V directly across the galvanometer circuit neces- 
sitated this, and the fusing system has been found to 
be entirely adequate. 


b. Decimal Point Section 


The decade voltage divider connected across the 
high-voltage terminals divides that voltage by integral 
powers of ten. The switch used for this purpose is 
the fourth switch to the right in the “V =” section on 
the main control panel (fig. 3). This switch also enters 
the appropriate power of ten in the Decimal-Point 
computer and lights the appropriate decimal point 
light on the main control panel (fig. 3). The rela- 
tionship between the decimal point section of figure 5 
and the “V=” dials of figure 3 is as follows. The 
10-' position causes the decimal point light preceding 
the first significant figure to light; hence the reading 


os % - J 


DECIMAL POINT SECTION 


The set-volts potentiometer. 


of the V=section, as shown in the figure, would be 
150. The 10° position causes the light between the 
first and second significant figures to light: hence 
the reading would be 1.50, etc. When the decimal 
point switch is rotated to the 10* position, a zero ap- 
pears in the fourth window of the V =section and the 
decimal point light to the right of that window will 
light. With the switch in this position the reading 
would be 1500. If a potential of 1500 V was now im- 
pressed across the decimal-point-section divider, the 
voltage to the right of G; would be 


1500.(10%) 


Tor = 0-150 volt, 


(6) 


which is equal to the voltage to the left, and the gal- 
vanometer would be at null. 


The possible errors introduced by the circuits in the 
set-volts potentiometer are discussed in appendix A. 


7.3. Computing Potentiometer 


The computing potentiometer (fig. 6) is basically a 
modified Lindeck-Rothe circuit connected to a voltage 


divider. The ammeter normally used as an indicator 
has been replaced by a 10 kf) resistor. Hence the 
output from the potentiometer is a voltage rather than 
a meter deflection. As is explained below, this voltage 
is operated on by the voltage divider comprising R 
and R: such that 


— Re yy. 
Eo=R M Vd, (7) 
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FIGURE 6. Computing potentiometer. 


where vq is the voltage to be measured by the poten- 
tiometer, and M is a multiplier determined by the 
values of Raq and Rp. 

The three decades shown under “V =” in figure 6 
are ganged to the corresponding decades in the set- 
volts potentiometer. It will be noted that when these 
decades are adjusted to the minimum value (V = 100) 
there is still 10 kQ of resistance in the circuit. It is 
this feature that allows R» of eq (5) to vary independ- 
ently of R:. . 

When both the “V =” and “‘S,,”” decades are adjusted 


to minimum, the circuit is as shown in figure 7a. If 


+h) 


the “V =” dials are adjusted to 325, for example, the 
resistance of R, becomes 32.5 k© and the result is 
shown in figure 7b. The value of Eo in figure 7b is 
equal to that of figure 7a divided by 325 (ignoring the 
decimal point). Thus the “V =” value is entered as a 
divisor. If the S, decade is adjusted to 4, the circuit 
is as shown in figure 7c, and Eo is four times that of 
figure 7b. Thus the value of S, is entered as a multi- 
plier, and the circuit of figure 6 meets the requirement 
of eq (5). 

The errors associated with the computing poten- 
tiometer are discussed in appendix B. 


7.4. Decimal Point Computer 


The decimal point is placed by means of the circuit 
shown in figure 8. Placement is determined by three 
switches, as shown. It will be noted that as the value 
of the Lindeck shunt selected is increased, the decimal 
point moves to the right. This is true also when the 
S, decimal point switch is changed to select a larger 
value. If the decimal point switch on the set-volts 
potentiometer is adjusted to a larger value, however, 
the decimal point on the digital display will move to 
the left. Thus the Lindeck shunt and S, decimal 
point switches serve to multiply by powers of 10, and 
the set-volts decimal point switch divides by powers 
of 10. In the actual circuit, provision is included for 
bypassing the S, switch when v; and v, are being 
measured. 


7.5. Switching Circuit 


The measuring circuit is connected to the voltages 
V1, V2, v3, and vs by means of double-pole, single-throw 
relays, as shown in figure 9. It will be noted that the 
movable and stationary contacts of any relay are sep- 
arately guarded. This is necessary to minimize the 
effects of leakage currents flowing between the con- 
tacts when the relay is deenergized. With the circuit 
as shown, any leakage currents of this type flow directly 
from the high-voltage supply and hence have no effect 
on the ratio. To implement the guarding shown in 
figure 9, it was necessary to modify commercially 
available relays to the form shown in figure 10. 


Roe lka Eo*306 


FIGURE 7. 


E 
2 


Computer circuit showing ““V="" and “‘S,”” adjustments. 


A: V=100, S,=1 
B: V=325, S,=1 
C: V=325, S,=4. 
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FIGURE 10. Switching circuit relay. 


7.6. Guard Circuit 


In keeping with personnel safety requirements, the 
cabinets of the console are connected to a low-im- 
pedance ground. Grounding the cabinets also pro- 
vides the measuring circuit with shielding. The 
measuring circuit guard is grounded also as shown 
in figure 9. The power supplies are isolated; one side 
of the high-voltage line is grounded (see fig. 9). When 
V1, V2, Or v3 is being measured, the side connected to 
v; is grounded. When v, is measured, the side con- 
nected to vs is grounded. 


8. Conclusion 


The incorporation of this console into the measure- 
ment facility at the Boulder Laboratories has resulted 
in a significant reduction in the limits of uncertainty 
associated with the calibration of volt boxes. The 
results of measurements made on a box of the type 
described in RP1419 are well within the predicted 
worst-case uncertaintly of +8 ppm (see appendix C). 


TABLE 1. Comparison of measurements made by the 
self-check method with those made on the console 


| 
Corrections in ppm Comparison 
|} in ppm 
— 4 
| | 
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| | 
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These results are shown in table 1. Column B shows 
the results of measurements made by the “self check- 
ing’ technique described in appendix II of RP1419. 
Columns C and D are the results of two independent 
measurements made by comparing the unknown 
(previously self-checked) volt box with the standard 
by means of the console and, in fact, reflect not only 
the console uncertainties, but the uncertainties in 
the correction to the standard ratio and any uncer- 
tainties introduced by the operator also. The com- 
parisons were made at 20 percent of rated voltage on 
each range to approximate the conditions under which 
the self check was made. Column E is the average of 
columns C and D. These conclusions are corrob- 
orated by the results of a second experiment that 
involved measurements of three master volt boxes 
in the following manner. 
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First V.B. 1 was measured against V.B. 2, 
then V.B. 3 was measured against V.B. 1, and 
finally V.B. 3 was measured against V.B. 2. 


In the first and third measurements, V.B. 2 was con- 
sidered the absolutely perfect standard. Inthe second 
measurement, V.B. 1 was considered the absolutely 
perfect standard. 

If the corrections, in ppm, to V. B. 1 resulting from 
the first measurement are called A; those to V. B. 3 
resulting from the second measurement are called B; 
and those to V. B. 3 resulting from the third measure- 
ment are called C, the following relationship should 
exist [5]. 

A+B-—-C=0, (8) 
provided the corrections to the three volt boxes did 
not change between measurements, and provided no 
errors are introduced by the measuring system—in 
this case, the console. To insure unchanging cor- 
rections in the volt boxes, measurements were made 
at reduced voltage. When the data from the measure- 
ments were applied to eq (8), the results were as 
shown in figure 14. 


The author is indebted to Joseph J. Barth, who de- 
ciphered the many schematics and wiring diagrams 
and reduced them to operating hardware. 
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10. Appendix A. Evaluation of the Set-Volts 
Potentiometer 


10.1. Kelvin-Varley Circuit (See Fig. 5) 


The rigorous mathematical analysis of a Kelvin- 
Varley divider is quite complex, and has appeared 
in the literature.[4] 

As the divider in the Set-Volts Potentiometer is to 
be used without applying corrections, the maximum 
effective uncertainty that will be introduced by the 
divider is the only thing that must be determined. 
Individual errors in the several sections are important 
only in their contribution to the overall error. As a 
recently calibrated digital voltmeter having combined 
uncertainties within 0.01 percent was available for 
the purpose, it was simpler and easier to measure this 
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settings. 
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FIGURE 12. Deviations in measured value of ratio caused by dis- 
agreement between value of applied voltage and value to which 
set-V dials are set. 


overall error than to calculate it. Accordingly, the 
900 settings on the set-volts divider were measured 
against the digital voltmeter, and the data were 
checked at 19 points with a commercial potentiometer 
(in which the combined uncertainties were within 
0.01%). The results of these measurements are shown 
in figure 11. As can be seen, the largest error intro- 
duced by the nonlinearity of the divider is roughly 
equal to one digit on the third (least significant) dial. 
When the applied voltage is adjusted to a value equal 
to one division in the third place higher or lower than 
the value to which the Set-Volt Potentiometer dials 
are set, the resulting difference in measured ratio, 
from that determined when the applied voltage and 
dial setting are equal, is smaller than the random 
errors encountered in the measurement. The results 
of such an experiment are shown in figure 12. 
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FIGURE 13. Decimal divider in set-volts potentiometer. 


10.2. Decimal Divider 


The decimal divider is shown in figure 13, which 
corresponds to the circuit of the Decimal Point Section 
in figure 5. In figure 13, the relationship of Eo to Ej 
for any given switch setting, n, is 





n-(——* 


ey ae ge (Al) 


where Ey and R,, Ro, . . . Rn are che nominal values. 
The true value, Ey, however, is ° 


Rx +r)E; 


> Rx(1 + rx) 
1 


Ey= Eo + €a)= 


Dividing eq (A2) by eq (Al) gives 


(1+) S$ Re 
1 


n be 


SRA + rx) 
1 


l+eqg= 


+r) ¥ Re-S Rl +70) 
= 2 2 
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= ’ (A4) 
> Rl + rx) 
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FIGURE 14. 


which reduces to 


(A5) 


Equation (A5) will be maximum when n=5 and re, rz, 
r4, and rs are opposite in sign to r;. In this worst case, 
assuming the corrections to all the resistors are maxi- 
mum, i.e., 0.0005,” eg is found to be 


€a = 0.001. (A6) 


11. Appendix B. Evaluation of Errors Asso- 
ciated With the Computing Potentiometer 


As explained in the text, the milliameter of the 
Lindeck-Rothe potentiometer has been replaced by a 
resistor (Ry) of 10,000 2 nominal value (see fig. 6). 
The voltage across R» will (nominally) be 10*, 10°, 
10°, or 107 times vg depending on which value of Ra 
is switched into the circuit. If this value, Vp,, is ap- 
plied across the divider circuit, R,, the nominal value 
of voltage, Eo, applied to the digital voltmeter would 
be found by the relation: 


(Bl) 


? The resistors used were individually measured and found to have corrections well 
within the limits of + 0.0005. 














CONSOLE-CAUSED ERRORS (A+B-C) IN ppm OF RATIO 


1 L 4 


30 45 60 75 150 300 450 600 750 
VOLTAGE RANGE IN VOLTS 





Console-caused errors as determined by intercompari- 
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where va is any one of the four difference voltages 
(v;. vo, v3, or vs), R, and Ry» are the displayed (i.e., 
nominal) values proportional to V and S, respectively, 
see eq (1), Ry is nominally 10,000 0, and R, is nom- 
inally either 1, 0.1, 0.01, or 0.001 0 as determined 
by the position of the Lindeck Shunt switch. How- 
ever, it is not practicable to measure Vr, without chang- 
ing its effective value (because of the loading effect 
of the measuring system) and hence disturbing the 
ratio of multiplication, R,/Ra. Therefore, to prevent 
this disturbance, the voltage is measured across 
R,+R,. With a circuit of this configuration, the 
voltage applied across R, is not Vx,, but Vr, +Ve,, 


and the actual value, Eo, is related to vq by the equation 


R; R,+Ri 
. 


Eo= (B2) 


In eq (B2) the primes indicate actual (as opposed 
to nominal) values. 

As providing automatic correction for the effect of 
R, in the numerator of eq (B2) would have appreciably 
increased the cost of the console, it was decided to 
treat the circuit as though the voltage applied across 
R, is Ve, even though in fact this is not the case. 
This will, of course, introduce an error in addition to 
that caused by the errors in R,, Ro, Ra, and Rp. It is 
the purpose of this appendix to determine the relative 
error (eo) in Eo, for any given va, resulting both from 
the errors in the component values and from the pre- 
sumption that eq (Bl) is the circuit equation. 

If the circuit of figure 6 is examined rigorously it 
will be found that the actual value of each of the com- 
ponents differs, although slightly, from the nominal 
value. If we consider each actual value (designated 
by a prime) to be comprised of the nominal value and 
the small difference between nominal and actual (i.e., 
Ri=R1+ra)), the circuit equation for determining 
the actual output voltage for a given vq is 


R.(1 + re) he + r)+ Ral sd , 


Eo= El +e)=F 4,7) Ra(1+ ra) 





(B3) 


where eo, ri, f2, Ta, and rp are the relative differences 
between the actual values, £6, Ri, Rs, Ri, and R; and 
the nominal values Eo, Ri, Ro, Ra, and Rp, respectively. 

Equation (B3) may be more readily analyzed if it 
is considered to be of the form 


E,(1 + eo) = T,0 + t:)T2(1 + te)va. (B4) 


If, in eq (B4), the following relations hold, 


R2(1 + re) 


Ti1+t)= Ril +r)’ 


(B5) 


Ro + To) + R.(1 + ra) 
R+ ra) 


T2(1 + te) = (B6) 





_R: 
Ry 
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aia 

then eq (B4) may be rewritten as 


T; 
and 


El te Ro qd a t:)(1 + te)vq. 


R, R, (B9) 


From a comparison of eq (B9) and eq (B1) it is imme- 
diately apparent that 


l+ep=(1+t)(1 + te). (B10) 


It is now possible to examine each of the component 
errors, ¢; and tz, individually; this possibility greatly 
facilitates the analysis. 

Let us first determine the relative error, te. 
divide eq (B6) by eq (B8), the result is 


If we 


1+ n= Ze 
b 


L+ro. 


1+Tre (B11) 


If |ra| and |r| are small (i.e., < 0.001), eq (B11) reduces 
approximately ° (i.e., to first degree terms) to 


_R 
te =—+1n—TIa.- 


“—_ (B12) 


In the circuit, the resistors R, and R, have been chosen 


such that |r| < 0.0004 and |r| < 0.0001. From figure 
6 it will be seen that 


Ra 


R, = 0.0001. (B13) 


Thus te will be within the limits 


\te|< (3) +\ral+|ro|= 0.0006. 
max 


Rp . (B14) 


The evaluation of t; may be handled in a similar man- 
ner. The division of eq (B5) by eq (B7) gives 


1+ re 


di eh or 
1 


(B15) 


3 If |ra| and |r| are each less than 0.001, the difference between the right-hand side of 
eq (B12) and the true value, 


—Ra 


1+nrm 
re. 


+ 
1+ra 


=i, (Fl) 


Tal» —Ta)| — 9 x 19-8 
l+ra 


and hence is so small as to not significantly affect the value of t2. 
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If |ri| and |r2|< 0.001, eq (B15) becomes 


t} =re—n. 


(B16) 


From figure 6, R; is seen to comprise many resistors 
in series. These resistors and switches have been 
chosen so that the actual value of the resistance added 
by any decade, on any position, does not differ from 
the nominal value by more than + 0.03 percent. Also 
the two 1 kQ resistors have been similarly chosen. 
Hence the circuit of R,; may be realistically repre- 
sented by eq (B17). (As R; is the total resistance of 
the network, it is represented by R; in the following 
analysis.) 


R(l+r)= > Rl + rv), (B17) 
k=1 


where R,, Ro, . . ., Ry are the resistors in the S, 
switch, the 1 kQ resistors connected to either end of 
that switch, and all other resistors that have been 
switched into the circuit by the settings of the three 
_ * dials. Equation (B17) may be rewritten in the 
orm 


Rl +r)=S Ret Reve. (B18) 
1 1 


(B19) 


R= Re, 
1 


eq (B18) may be divided by R; to give the equation 


r= BS Rin=(S Rare) /(SRe)- 20) 
es 1 1 


That is, 7 is dependent upon the sum of the differences, 
in ohms, between the actual and nominal values of 
the resistors Ry, in ratio with the sum of the nominal 
values. If, as previously stated, all values of |r;| are 
within some limit, r, then the value, r, as substituted 
for |r;| may be factored out of eq (B20) to yield the 
bounds to 7; i.e., 


DR 


1 


In| <r- =r. (B21) 


5 Ri 


Therefore, the bounds to 7; (i.e., to the relative differ- 
ence between the actual and nominal values of R;, 
in fig. 6) are 0.0003. Similar reasoning* can be used 


to determine the bounds to the relative error of R» 
(fig. 6) with the result: 


\r2| < 0.0003. (B22) 


Hence, the bounds to ¢; are found, using eq (B16), to be 


\t1] < 0.0006. (B23) 
However, the validity of this value is questionable 
since the interdependence of the errors of R; and R2 
was not considered in the analysis. A more rigorous 
analysis results when the resistance network R; is 


considered to comprise the independent resistances 
Rs and R» such that 


R3+ R2=R,. (B24) 


If this is done, eq (B5) becomes 


RA1 + re) 


ia ees ea 





(B25) 


The determination of the error introduced by R3 and 


_ Ry may be found by operating on the relation 


R2 
R3 + Re 


R2(1 + re) 


qd laa + r3)+ R2(1 +r) 





(B26) 


Dividing eq (B26) by R2/(R3 + Re) gives 


(R3 + R2)(1 + 12) 


to 
er R31 + rs) + Rol +72) 





(B27) 


If |rs| and |r2|< 0.001 (this has already been demon- 
strated as the reasoning used in the determination of 
r applies with equal validity to rs), eq (B27) may be 
approximated by (the reasoning of footnote 4 applies) 


R3+ Re 
R2.+ R3(1 + rz — r2) 





1+e,= (B28) 


being more than adequately accurate for the purpose 
intended. Subtracting 1 from both sides allows eq 
(B28) to be reduced to the more usable form, 


a R3(r2— rs) 
R.+ R3(1 + r3— re) 





er 


(B29) 


From eq (B29), |e,| approaches the value 


F2—_ U3 
l+r3—f 


* As Rz is in parallel with the input impedance of the digital voltmeter, the actual re- 
sistance may differ from nominal by as much as 0.09 percent unless compensating tech- 
niques are used. These techniques, in the form of trimmer resistors, were used on Ro. 
Hence the resistances shown are the effective values including the input impedance of the 
digital voltmeter, adjusted to within + 0.03 percent (300 ppm). 
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as a maximum limit as the ratio R3/R2 increases without 
bound. Using bounds of +0.0003 for r3 and rp 


0.0003 + 0.0003 _ 
‘1— .0003 — .0003 





ler| < 


(B31) 


le,| < 0.0006. (B32) 
A comparison of statement (B32) with statement (B23) 
shows that (B23) is indeed valid. 
Referring again to eq (B4), the value, eo, is quite 
approximately determined by the relation 
eo = ti tte (B33) 
provided |t,| and |t2| are < 0.001. That these values 
are less than 0.001 has been demonstrated. There- 
fore e9 may be determined from eq (B33) (using the 
limits of t; and tz listed in statements (B14) and (B23)) 
to be 
|eo|<|0.0006|+|0.0006| = 0.0012. (B34) 
This figure, found by direct addition, is certainly the 
pessimistic limit in the worst case (i.e., where 
R,/R2— ~~). A more realistic appraisal might be to 
combine the different independent errors by the square 
root of the sum of the squares, as they are as likely 
to compensate for one another as to add. If this is 
done using eq (B12). and (B16), then with high 
probability 





(B35) 


|eo|< Vrit+ r3+(Ra/Ro)? + m+ ra 


the result is 


|eo|< 0.0006. (B36) 
The actual limits. of eo will lie somewhere between 
these two extremes. 

Hence the premise of eq (B1) is found to be uncer- 
tain by about+ 0.1 percent of Eo. 

Total Measurement Circuit Uncertainty. An ad- 
ditional uncertainty of +0.2 percent (em) of the 
measured voltage is introduced by the digital volt- 
meter. Also the digital display has an inherent un- 
certainty of tone digit in the last place. Because 
of the way in which the display is used, this + one 
digit causes an absolute uncertainty (eq) of +1 ppm of 
the ratio (of the unknown volt box) in the worst case 
for each displayed value. 


12. Appendix C. Evaluation of the Measure- 
ment Uncertainty of the Volt Box Console 


This evaluation does not include the effect of the 
uncertainties associated with the corrections to the 


Standard Volt Box. 


The equation of correction for the unknown volt 
box as used by the operator of the console is 


X {1 + (x)10-°] = Xi{1 + (us + Vi + V2 + V3.4 V4)10-4, 
(C1) 


where X, is the nominal ratio of the unknown, x is 
the correction to the unknown in ppm, V;, V2, V3, and 
V, are the computed results in ppm, and yp; is the cor- 
rection to the standard volt box in ppm. 

The error in V,, V2, V3, or Vs is closely approxi- 
mated by the relation: 

e ~ ey— ey, (C2) 
where e, is the error resulting from the computer 
and ey is the error caused by the set-volts poten- 
tiometer. The difference in errors shown in eq (C2) 
is the reason that the zener supply for the set-volts 
potentiometer need only be calibrated against the 
digital voltmeter. If the digital voltmeter is-in error 
(absolutely) by some value, the set-volts potentiom- 
eter will be in error by the same value, and the two 
errors will subtract to zero. 

The uncertainties in e resulting from the dividers 
in the set-volts potentiometer, the resistors in the 
computing circuits, etc., must be considered as the 
signs of these uncertainties are not known. 

e, and ey are the actual values of the errors and, 
of course, are not known. It is possible to determine 
the approximate bounds within which these values 
must fall. 

The error, e,, is comprised of the contributing er- 
rors éo (the error in the voltage measured by the digital 
voltmeter), @m (the error introduced by the digital volt- 
meter circuits), and eg (the uncertainty in the last 
digit of the display) (see appendix B). These com- 
ponents have been determined to have the limits 


|eo|< 0.0012 of voltage measured by the digital 
voltmeter 


lem|< 0.0002 of voltage measured by the digital 
voltmeter 


ea==+1 ppm of measured ratio. 


The values of the limits of e9 and e» are functions 
of the voltage that is measured by the digital voltmeter 
and can be combined. Since eo and e» are independ- 
ent errors, their uncertainties can be combined by 
the relation: 


le, |< Vej+ e3, (C3) 


to give the highly probable limits to this portion of 
ey. Thus 





le}|< V(0.0012)? + (0.0002)? (C4) 


\e’|< 0.0012. (C5) 
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However, there is a possibility that the component 
errors could, under some circumstances, be of such 
sign that they would add, giving a value approximating 
the simple sum of the limits of uncertainty. Under 
these conditions 


le*|< 0.0014. (C6) 

The value of eq must be treated separately. 

The error ey, is also comprised of component 
errors. These are ea, the error resulting from the 
divider of figure 13, (see appendix A) and the errors 
introduced by the Kelvin-Varley divider. 

The value of e, has been determined to be within 
the limits +0.0001 of the applied voltage, V. The 
errors resulting from the Kelvin-Varley divider are 
somewhat obscure, but conservatively may be as- 
sumed to lie within the limits +2 ppm of the measured 
ratio. 


The accuracies of the values of V;, V2, V3, and Vs_ 


are directly dependent upon the accuracy with which 


the applied voltage can be measured (see eqs (2) and 
(3) in the text). However, the error in the measure- 
ment of the applied voltage will be the same for each 
value, and, though the sign is not known, it will be the 
same also. Hence the uncertainty of the measured 
ratio will very probably lie within the limits 


|u| Ss [Vier + Vrer+ Ve. + V se’) 


+(V;tV.t+V3+ V sea t+ 4eat exr| (C7) 
where ex, is the value of the limits of the error intro- 
duced by the Kelvin-Varley divider in ppm of ratio. 
When the limits of the different component uncertain- 
ties are substituted in eg (C7), the limits of u are found 
to be approximately 


|u| <8 ppm. (C8) 


(Paper 70C3-229) 
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Specimens of ten binary alloys of nickel and iron containing 3, 5, 10, 16, 20, 30, 36, 50, 57, and 81 
percent Ni and of the metals nickel and iron were immersed for 7 months in 265 gallons of city water 
to which had been added 3 percent by weight of sodium chloride. 

Cathodic and anodic polarization curves of the specimens were obtained at about 22 periodic 
intervals throughout the exposure period. Corrosion currents, calculated from currents at breaks in 
the curves, were converted to corrosion rates expressed as weight losses by applying Faraday’s law. 
The calculated weight losses were in reasonable agreement with the actual weight losses. 

The polarization rates, AV/AI, of the specimens calculated from cathodic polarization curves 
(cathodic control prevailed) at several intervals throughout the exposure period were plotted on log- 
arithmic coordinates with respect to the corrosion current densities calculated from breaks in the same 


curves. 


slope of —1. 


Similarly, the averages of these AV/A/ values were also plotted on logarithmic coordinates 
versus the actual corrosion rates expressed as weight losses. 


Both plots conformed to the theoretical 


Key Words: Corrosion rates, Ni-Fe alloys, polarization techniques, polarization rate, polarization 
resistance, polarization circuits, Ni-Fe pitting. 


1. Introduction 


A few years ago, it was demonstrated experimentally 
that the corrosion rates of a series of ferrous alloys, 
containing chromium additions up to 18 percent and a 
similar series with 3 percent of silicon in addition to 
the chromium, exposed to a 3 percent sodium chloride 
solution, could be calculated from currents associ- 
ated with changes of slope (breaks) in polarization 
curves [1].!. The sensitivity of the method was shown 
by the fact that the weight loss of the 18 Cr alloy 
attributed to corrosion was measured satisfactorily 
even though the apparent area affected by corrosion 
was less than 1 percent of the entire surface. The 
polarization curves of the alloys showed that the known 
critical change in corrosion rate (expressed as weight 
loss) for this range of chromium occurred when the 
corrosion reaction changed from cathodic to anodic 
control, between 11.4 and 14.3 percent Cr. 

The reason for the laboratory measurements de- 
scribed in this paper, using binary alloys of nickel and 
iron, was to explore further the practicability of polari- 
zation techniques as methods of evaluating rates of 
corrosion. The corrosion rates of several binary alloys 
between 3 and 81 percent nickel, including the un- 
alloyed metals, were measured by the above technique 
and also by a method which has been called “‘polariza- 
tion resistance” but is referred to by this writer as 
polarization rate AV/Al. 


'Figures in brackets indicate the literature references at the end of this paper. 


The significance of the slope, AV/AI, of the control- 
ling polarization curve was observed by the author and 
related to the break method [2] but was never further 
pursued until after Skold and Larson [3] found the 
method favorable. They observed that a linear rela- 
tionship existed when corrosion rate (expressed as 
weight loss) was plotted versus AV/AI on logarithmic 
coordinates. Shortly thereafter, Stern provided a 
theoretical basis for the polarization resistance method 
[4]. These methods among others have also been 
discussed by Phelps [5]. The writer has found it 
desirable to use both techniques simultaneously as an 
aid in interpreting and checking data [6, 7] and this 
procedure was followed in the measurements about 


to be described. 


2. Experimental Procedure 


2.1. Preparation and Exposure of Specimens 


The binary alloys were forged and heat treated by 
the Armco Steel Corporation. Their compositions 
according to analysis are shown in table 1. Forgings, 
0.875 in. X 0.875 (2.2 cm X2.2 cm) in cross section, 
were machined into specimens 0.625 in. (1.59 cm) diam 
X 12 in. (30 cm) long. A hole, 0.25 in. (6.3 mm) diam 
centered 0.32 in. (8.1 mm) from one end, was drilled 
through the specimen to permit insertion of a Nichrome 
wire, No. 30 B&S gage, 0.010 in. (0.25 mm) diam, which 
served as a loop to hold the specimen in suspension 
below the surface of the corrodent. All edges were 
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rounded, particularly at the hole, in order to assure 
accuracy in weight loss measurements. This precau- 
tion also prevented nicking of the suspension wire and 
provided good electrical contact as this wire also con- 
ducted the polarizing current. 

The method used in suspending the specimens and 
the reasons therefor were previously described [1]. 
The Nichrome was always cathodic to the specimens, 
even to the nickel specimen, and galvanic action, if any, 
was insignificant as previously observed with the Fe-Cr 
alloys. 

Before exposure, the specimens were degreased, 
scrubbed with scouring powder, rinsed, rubbed with 
fine-grit abrasive, scrubbed under running hot water, 
dried by air-blast and then weighed to the nearest 
2 mg. 

Twelve specimens, each of different composition 
(table 1) were exposed to Washington, D.C. city water 
to which was added 3 percent by weight of sodium 
chloride. The salt water was contained by an open 
cylindrical vat, 66 in. (16.7 dm) diam. The depth of 
the water was maintained at 18 in. and the volume of 
water about 265 gal. The specimens, about 13 in. 
(33 cm) apart, were staggered and suspended in two 
rows across the center of the vat, the tops of the speci- 
mens being 2 in. (5 cm) below the water surface. 
Electrical contact with a given specimen was made 
through a terminal strip mounted on the outside of the 
vat from which wires were permanently connected to 
the Nichrome suspension loops by battery clips. The 
temperature of the water was not controlled but was 
measured regularly during exposure. 


TABLE 1. Composition of alloys 








Compositon, percent* 





Specimen 


Iron 


Nickel 


100 
97 
95 
91 
85 
81 


| 
| 
| 
COTA | 
| 
| 
| 











*Analyzed by x-ray fluorescence techniques for iron and nickel in each specimen. Iron 
and nickel were determined separately. Values as measured were within +0.5 percent 
and are rounded off to the nearest 1 percent. Two specimens, | and 11, were analyzed 
bya — spectrochemical method for impurities, resulting in the following: Specimen 

‘1 had Ni, Cu, Si, Ti, each between 0.01—0.1 percent, Al< 0.01 percent. Specimen 11 
had Si mane 0.1-1.0 percent, Cu, Co, Ti, each between 0.01—0.1 percent, Al, Mn each 
< 0.01 percent. 


2.2. Electrical Measurements 


Potentials were measured with reference to a satu- 
rated calomel half-cell. Contact between the salt 
water and the half-cell was by an agar-salt bridge in the 
form of a 5-ft (152 cm) length of flexible tubing contain- 
ing a mixture of agar and potassium chloride. The 
electrical continuity of the agar-salt bridge was assured 


‘vat. 


by a cotton cord saturated with potassium chloride. 
The cord extended the entire length of the tube 
enveloped by the agar-salt mixture. The end of the 
tube was placed into the water through a sleeve 
mounted about 6 in. (15 cm) from the wall inside of the 
The distance between the tip of the salt-bridge 
and a specimen varied between 18 in. (46 cm) and 40 in. 
(101 cm), depending on the location of the specimen. 
When electrical measurements were not in progress, 
the salt-bridge was removed from the water. When 
polarization measurements were being made, aux- 
iliary electrodes consisting of two strips of zinc sheet, 
4 in. (10 cm) wide X 18 in. (46 cm) long, were hung over 
the wall of the vat. 

Measurements were made on about 22 occasions 
during the exposure period of 210 days. The measure- 
ments involved corrosion potential and the running 
of cathodic and anodic polarization curves. For a 
given specimen, an interval of several hours, some- 
times 24 hr, was permitted between cathodic and 
anodic polarization runs, not necessarily in that order, 
but usually so. 

Corrosion potentials were measured with a high 
resistance voltmeter and polarization was measured 
and recorded with a two-pen (current and potential) 
strip-chart recorder. The length of time required for 
obtaining a polarization curve varied from about 15 
to 30 min. Most of the polarization data were obtained 
by using the galvanodynamic method and about 25 
percent using a potentiodynamic method for compari- 
son. The schematic circuit diagrams for these 
methods are shown in figures 1 and 2, respectively. 
Both circuits incorporate the Holler bridge for balanc- 
ing out undesirable voltage drop (IR) from the recorded 
potentials [8]. In figure 1, the polarizing current is 
steadily increased by varying the applied voltage at a 
linear rate by means of the ten-turn motor-driven poten- 
tiometer R,; with resistor R3; set for a predetermined 
range of polarizing current. In figure 2, initially the 
applied voltage between the specimen and the refer- 
ence electrode is manually adjusted by setting R4 to a 
value equal to the corrosion potential after which R4 
becomes motor-driven and sweeps the range of poten- 
tial required for polarization while Rs; is set at a com- 
paratively low value of resistance so as to not unduly 
limit the applied current. In viewing the recorder 
charts, the break in the curve is more evident on the 
potential plot than on the current plot with the galvano- 
dynamic method, while the reverse is true when the 
potential sweep circuit (fig. 2) is used. 


4 


2.3. Removal of Corrosion Products 


Following removal from the salt water, the specimens 
were held under running hot water and scrubbed with a 
stiff fiber-bristle brush. They were then placed into 
an ultrasonic cleaning tank (for 10 min) containing a 
10 percent solution, at 150 °F, of ammonium citrate 
made alkaline (pH 9-10) by the addition of ammo- 
nium hydroxide. Then, the specimens were again 
brushed under running hot water, dried under an air 
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FiGuRE 1. Polarization circuit-galvanodynamic method. 
R, = 250 — 10 turn (motor dfiven—0.2 RPM), R2= 10010 turn, Rs, 10,0002, X = 1000 
— 10 turn, @= D= 100,0002. 
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FIGURE 2. Polarization circuit —potentiodynamic method. 
R,=500 — 10 turn (motor driven—variable speed—reversible), R2=100Q 10 turn, Rs 
= 10,0008, Ry=25Q—10 turn (motor driven —0.2 RPM reversible), X = 1000 — 10 turn, 
Q=D= 100,0000, B= motor (controlled by voltage differential). 


_blast and weighed to the nearest 2 mg. After repeat- 
ing this procedure the third time, there was no signifi- 
cant change in weight; the final weights were then 
subtracted from the original weights in computing the 
actual metal loss attributed to corrosion. 


3. Results and Discussion 


On the 210th day of exposure the specimens were 
removed from the salt water. Although the object of 
the experimental work was primarily that of evaluating 
the corrosion rate measuring techniques, the data in 
table 2 show how the results to be described are related 
to the apparent areas corroded and to the depths of the 


pits. Note, that all apparent corrosion on the nickel 
(specimen 12) seems to have been confined to one pit. 

Polarization curves from data on recorder charts 
obtained on the 202d day of exposure are shown in 
figures 3 and 4 on semilogarithmic and rectangular 
coordinates, respectively. While such data, ordi- 
narily, are shown on semilogarithmic coordinates, the 
rectangular plot helps to verify the break in the curve 
and the value of applied current at which it occurs. 
The absence of a break on the latter would place doubt 
on the extrapolated semilogarithmic value. The 
values of J, (break in the cathodic curve) and of Iq 
(break in the anodic curve) chosen in calculating the 
corrosion currents (table 3) are marked /, and J, in 
figures 3 and 4 and are indicated by vertical dashes 
(where dashed lines intersect). It will be noted that 
there is reasonably good agreement between figures 
3 and 4 for a given specimen in the values of /,, but 
apparently not as good agreement in the values of J). 


TABLE 2. Appearance of the binary alloys of nickel and iron after 
cleaning which followed 210 days of exposure to salt water 








Apparent 
area 
corroded 


Specimen Remarks, pit depths 





Percent 
1 (100 Fe) 85 General corrosion (0-3 mils*). 


2 (3 Ni) 70 General corrosion (0-3 mils), except for one pit (5 mils). 


3 (5 Ni) 80 General corrosion (0-3 mils). 


4 (10 Ni) 75 General corrosion (0-3 mils). 


5 (16 Ni) 65 General corrosion (0-5 mils). 


6 (20 Ni) 50 General corrosion (0-8 mils). 


7 (30 Ni) 40 Lower 4.5 in. (11 cm) (opposite suspension end) of specimen 
length not corroded, otherwise — general corrosion (0-8 mils). 


8 (36 Ni) Spotty corrosion to 10-mil depths. 


9 (50 Ni) 8 in. (20 cm) of the specimen measured from the lower end was 
not corroded. 

Spotty corrosion (8-12 mils) on the upper 4 in. (10 em). 
45 mils in depth at the top end. 


Pits to 


10 (57 Ni) Scattered pits (2-6 mils) in about 24 places. One pit (50 mils) 
at the top end. 


11 (81 Ni) Three pits (5-25 mils) on the surface. One pit (50 mils) inside 
of the mounting hole. 


12 (100 Ni) No corrosion, except for one pit (60 mils deep) inside the mount- 


ing hole. 











*mil=0.001 in. (25.4 y). 


The values of J,, pertaining to these data, selected from 
the rectangular coordinates are considered the more 
significant, except for specimens 11 (81 Ni) and 12 
(100 Ni). The actual values of J, and J, were obtained 
before compressing the scales (particularly the ordi- 
nates) as shown in figures 3 and 4. The data for all 
specimens are shown on one figure so that the reader 
can get an overall picture of what the curves portray for 
different rates of corrosion. When calculating the 
corrosion current (see the equation, footnote b, table 
3), the smaller current, whether it be J, or Jg, is the 
more significant (controlling) and fortunately can be 
extrapolated with the greater degree of accuracy. 
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FIGURE 3. Polarization curves of iron, nickel, and binary alloys 
of nickel and iron obtained on the 202d day of exposure to city 
water having added 3 percent by weight of sodium chloride. 

O anodic, @ cathodic. 

NOTE: Values of J, and J, actually used for the calculations (table 3) are marked /, and 
1, (figs. 3 and 4), for example, /, (fig. 3) is used with /, (fig. 4) for the 100 Fe specimen. In 
the case of the 100 Ni specimen, /, and /, of figure 3 are used, etc. The ordinates (figs. 
3 and 4) are compressed to save space. In going up the scale, the potentials become less 
negative. Thus, for the 30 Ni specimen, the corrosion potential is about — 0.58 V (see fig. 4). 


All the values of J, and J,, calculated corrosion cur- 
rents, calculated weight losses, and actual weight 
losses for all specimens are tabulated in table 3. In 
calculating weight losses, the electrochemical equiva- 
lent for ferrous iron was chosen for the alloys with less 
than 50 percent of nickel, and that of nickel for the 
alloys having 50 percent or more of nickel. Actually, 
the equivalent value chosen is not significant here as 
the difference (about 5%) in results obtained is within 
the limit of accuracy with which the polarization curves 
can be evaluated. 

In table 3, beginning with the 78th day of exposure 
and on through the 117th day, all data were obtained 
using the potentiodynamic method of polarization. It 
will be noted that the data pattern fits in quite well with 
the pattern before and after this period obtained using 
the galvanodynamic polarization technique. For 
specimens | through 6, the change from J,=10 I, to 
lesser values of J, is not believed to be significant. 
The potentiodynamic method offered no advantage in 
obtaining these data and required more instrumentation. 


V-S.GE. 





POTENTIAL, 











QO .02 .04 .06 .08 
APPLIED CURRENT, mA 


FIGURE 4. Same data as in figure 3, plotted on rectangular 


coordinates. 
O anodic, ® cathodic. See “Note” figure 3. 


Figure 5 shows how the corrosion rates vary with the 
nickel content and how the calculated and actual 
values compare. The biggest reduction in corrosion 
rate occurred in changing from 50 Nito57 Ni. Forthe 
particular environment, nickel in amounts up to 20 
percent seem to offer no advantage from the stand- 
point of corrosion resistance. Pettibone [9] demon- 
strated an advantage in the corrosion resistance of 
36.53 Ni steel over mild steel in four sea water environ- 
ments for exposure periods varying from 5 to 15 years. 
It has also been reported that the corrosion rate of a 
26 percent Ni alloy of iron and nickel was about one- 
third that of wrought iron in sea water as well as in 
the atmosphere [10]. 

When specimen 8 (36 Ni) was initially exposed, the 
corrosion current was about one-eighth that of the 
specimens with lesser amounts of nickel (table 3). 
Between the 4th and 7th days, the corrosion potential 
of specimen 8 changed from —0.355 V to —0.540 V 
{potentials are not shown in the table) and the corro- 
sion current increased fivefold. On specimen 9 (50 
Ni), a low rate of corrosion prevailed through the 34th 
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TABLE 3. Corrosion weight losses calculated from polarization TABLE 3. Corrosion weight losses calculated from polarization 
data obtained on a series of alloys containing nickel and iron data obtained on a series of alloys containing nickel and iron 
exposed to a 3-percent sodium chloride solution for 7 months exposed to a 3-percent sodium chloride solution for 7 months— 
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TABLE 3. Corrosion weight losses calculated from polarization 
data obtained on a series of alloys containing nickel and iron 
exposed to a 3-percent sodium chloride solution for 7 months— 
Continued 


TABLE 3. Corrosion weight losses calculated from polarization 
data obtained on a series of alloys containing nickel and iron 
exposed to a 3-percent sodium chloride solution for 7 months— 
Continued 
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* See the text. 

> io=I1p-Iql (lp + IQ). 

© Calculated weight loss (g)=Kti, K=electrochemical equivalent (for specimens 1-8, 
assumed K=2.8938 x 10-* g/coulomb; for specimens 9-12, assumed K = 3.0409 x 10-4 
g/coulomb; i= average of corrosion currents (A) for the exposure period (t in sec) between 
successive readings). The values of io at the beginning and end of exposure are taken as 
the initial and final values, respectively, as calculated. 

* Assumed. 


NOTE: All speci were completely submerged and exposed at the same time in an 


open vat containing about 265 gal of salt water. Mean solution temp., 73 °F; min, 64 °F; 
max, 81°F. Area of specimen, 1/6 ft? (155 cm?). 
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FIGURE 5. Effect of nickel alloyed with iron on the corrosion rates 
of alloys exposed for 7 months to Washington, D.C. city water to 
which was added 3 percent by weight of sodium chloride. 

@ based on the actual weight loss for 210 days. 


O based on the cumulative weight loss as calculated from the breaks in polarization 
curves (table 3). 


day (table 3) after which the corrosion potential 
changed from — 0.343 V to — 0.495 V and the corrosion 
current increased sevenfold. The corrosion current 
pertaining to specimen 9 was minimum when its corro- 
sion potential was the least active (—0.285 V) on the 
7th day of exposure. For a comparison of the corro- 
sion potentials of all alloys on the 202d day of exposure, 
see figure 4 (at zero current). Perhaps, the still water 
surrounding the specimens, where a greater chance 
for oxygen depletion occurred, may have prevented a 
continuation of the film-forming tendencies initially 
exhibited by the 36 and 50 Nialloys. In moving water, 
the oxide films would probably have been more stable 
and might have prevailed. 

The relationship between the two corrosion rate 
measuring techniques is shown in figure 6. The line 
through the points is drawn with a slope of —1, the 
theoretical basis for which was given by Stern [4]. 
The values of AV/A/ are based on the change in poten- 
tial caused by the increase of applied current from 
zero to a value slightly larger than the current /,, as 
indicated by arrows in figure 3 when the corrosion 
current is zero or negligible as indicated by the in- 
crease polarization rate. Selecting values of applied 
current smaller than /, could result in erroneous values 
of AV because the changes in potential might be too 
small for reliable measurement. Once a relationship, 
as in figure 6, is established for a given environment, 
corrosion rates can be approximated simply by measur- 
ing AV/AI. 

As the corrosion currents were controlled by the 
cathodic reactions, almost without exception (table 3), 
the values AV/Al were calculated from the cathodic 
polarization curves. The relationship between the 
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FIGURE 6. Logarithmic relationship between the corrosion current 
density, calculated from breaks in cathodic and anodic polariza- 
tion curves, and the controlling polarization rate (cathodic) of 
nickel, iron, and binary alloys of nickel and iron. 


Sublegend: Seven values are shown for each specimen based on measurements made 
at periodic intervals during the 210-day exposure period. 
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C) Specimen 10, 57 Ni. 
A Specimen 11, 81 Ni. 


@ Specimen 12, 100 Ni. 


polarization rates AV/AI and the actual corrosion rates 
(weight losses) of the 12 specimens is shown in figure 7, 
the line through the points again being drawn with 
slope of — 1. 


4. Summary 


Ten binary alloys of nickel and iron (series of 3 to 
81% Ni) and the metals nickel and iron, 12 specimens 
in all, were exposed for 210 days to about 265 gals 
of still Washington, D.C., tap water (64-81 °F) to which 
was added 3 percent by weight of sodium chloride. 
Corrosion currents were calculated periodically from 
breaks in polarization curves, averaged for each of 
the 22 periods and converted to cumulative weight 
losses by Faraday’s law. The calculated weight 
losses differed from the actual weight losses (resulting 
from corrosion) by from 1 to 61 percent, averaging 22 
percent for the 12 specimens. The iron lost about 85 
times as much weight as did the nickel while the ratio 
of the respective areas apparently corroded was 
about 4000 to 1. Corrosion of the nickel seemed 
to be confined to one pit about 0.064 in. (1.6 mm) diam. 
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FIGURE 7. Logarithmic relationship between corrosion rate based 
on actual weight loss and the controlling polarization rate of 
nickel, iron, and binary alloys of nickel and iron. 

Sublegend: Each point pertains to a specimen and the polarization rate is the average of 
seven measurements made at periodic intervals during the 210-day exposure period. 


@ Specimens | through 7, 0-30 Ni. 


) Specimen 8, 36 Ni. 


D Specimen 9, 50 Ni. 
| Specimen 10, 57 Ni. 
\ Specimen 11, 81 Ni. 


@ Specimen 12, 100 Ni. 


Polarization rates, AV/A/J, of the 12 specimens, 


referred to by most investigators as polarization resist- 
ance, were plotted on logarithmic coordinates against 
corrosion current densities calculated from the same 
curves by the break method and found to fit a curve 
having the theoretical slope of —1. Similarly, the 
actual weight losses of the 12 specimens plotted on 
logarithmic coordinates versus their average polariza- 
tion rates also revealed a slope of —1. Thus, once 
having established the relationship between corrosion 
current density (from breaks in polarization curves) 
and polarization rates, AV/A/J, corrosion current 
densities can be approximated for other alloys in the 
particular environment by simply measuring AV/AI. 


The primary purpose of the laboratory work was to 
evaluate two polarization techniques as methods for 
measuring rates of corrosion:in an aqueous solution. 
Of secondary interest to the author was the observa- 
tion that the biggest benefit in corrosion resistance 
occurred in changing from 50 to 57 percent Ni. Both 
the calculated and the actual weight losses revealed 
a reduction in corrosion of over 80 percent for the 57 
Ni alloy over that of the 50 Ni alloy. 


The author gratefully acknowledges the cooperation 
of the personnel of the Research Laboratory, Armco 
Steel Corporation, Baltimore, Md.. who arranged for 
the fabrication and prepared the alloys. 
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Creep tests were made at 700, 900, and 1200 °F (644, 755, and 921 °K) on a nickel-4.2 percent 


aluminum alloy (Duranickel) initially as cold-drawn and as cold-drawn and age-hardened. 


Micro- 


structural changes occurring during creep are correlated with the flow, fracture, and ductility of the 


specimens. 


Light and electron microscope observations indicate that increases in creep strength 
to] I ba} 


were associated with the formation of rather straight discontinuous slip bands, a well defined subgrain 
structure, small precipitates of Ni;Al, and a general distribution of these precipitates throughout the 


grains. 


Creep strength at 1200 °F was relatively unaffected by prior cold-drawing or aging. 


Inter- 


crystalline cracking, accompanied by low ductility, was evident at 900 and 1200 °F and slow creep rates: 
whereas transcrystalline cracking and high ductility values were observed at 700 °F. 


Key Words: Cold-drawing, creep, nickel alloys, Duranickel, nickel-aluminum precipitates, high 


temperatures, deformation. 


1. Introduction 


A comprehensive investigation designed, to evaluate 
the influence of stress, strain, strain rate, temperature, 
and structure on the rheological properties of a precipi- 
tation-hardening nickel-4.2 percent aluminum alloy 
(Duranickel*) has been in progress for the last few 
years at the National Bureau of Standards. The 
tensile and creep properties of annealed and annealed 
and aged Duranickel have been previously published 
[1,2].!. The present phase of the investigation 
describes (1) the influence of cold-drawing on hardness 
and creep behavior; (2) the effect of prior aging on 
creep characteristics; (3) the influence of creep on 
aging; (4) the role of creep rate and microstructure on 
crack formation and fracture. In addition, data, 
previously published [2] describing the temperature 
dependence of creep behavior of annealed and an- 
nealed and aged Duranickel are included for compari- 
son purposes. 


2. Material and Testing Procedure 


The chemical composition of the metal used in this 
investigation is given in table 1. All the bars were 
processed from the same heat as that previously used 
[1, 2] and were supplied by the Huntington Alloy 
Products Division of the International Nickel Com- 
pany in the form of 7/s-in. round rods, approximately 
15-ft long. The bars were hot-rolled, descaled, and 


*Trade name International Nickel Company. 
' Figures in brackets indicate the literature references at the end of this paper. 


annealed at 1725 °F (1213 °K), then cold-drawn by the 
manufacturer to give a tensile strength of 127,000 psi 
and a hardness of 27 Rockwell “C’. Specimens, 
Yq in. in length, were cut from the bar and aged in 
vacuo for 16 hr at temperatures ranging from 700 to 
1250 °F (644 to 948 °K). The specimens age-hardened 
above 900 °F (755 °K) were cooled at a rate of 25 °F/hr 
to 900 °F, then furnace cooled to room temperature. 


TABLE 1. Chemical composition (percentage by weight) of the 


Duranickel used 





Determinations made at the National Bureau of Standards. 


Specimens age-hardened at 900 °F and below were 
also furnace cooled. Room-temperature hardness 
surveys indicated maximum hardness at 1100 °F (866 
°K) was attained with a 16 hr aging time (fig. 1C). 
Above this temperature at this time the metal appeared 
to be overaged (fig. 1A). In order to ascertain the 
effect or rate of cooling on the initial structure and 
hardness of the cold-drawn metal, specimens were 
heated in air for 16 hr at temperatures ranging from 
700 to 1200 °F (644 to 921 °K) and water quenched 
(fig. 1A). Room temperature hardness surveys in- 
dicated the hardness values for each temperature were 
lower than those obtained by furnace cooling. Another 
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FIGURE 1. Effect of aging temperature and time and rate of cooling on the hardness of Duranickel. 
FC —furnace cooled; WQ—water quenched. 


group of cold-drawn specimens was aged at 1100 °F 
for 16 hr and furnace cooled. Then they were heated 
in air to temperatures of 700, 900, 1100, and 1200 °F 
and held for 16 hr before water quenching. This 
procedure was used in order to obtain more stable 
microstructures in the specimens. Room-temperature 
hardness values indicated a softening effect for the 
cold-drawn and aged specimens (fig. 1B). The same 
general trend was also observed for the annealed and 
aged specimens (fig. 1B). 

Cold-drawn rods were cut into specimens 6%/4-in. 
in length for creep testing. These were divided into 
two groups: one group to be tested in the cold-drawn 
condition, the other group to be tested in the aged 
condition (aged at 1100 °F for 16 hr and furnace cooled). 
Room-temperature tensile properties were determined 
by the authors for the annealed, cold-drawn and aged 
specimens (table 2). 

The creep specimens had a reduced section of 0.252- 
in. diam over a 2-in. gage length and the ends were 


TABLE 2. Room temperature tensile properties of Nickel-4.2 percent 
aluminum alloy 











Yield 
strength 
(0.2% offset) 


Tensile 
strength 


Elon- 
gation 





Percent | Percent 
47 77 
20 48 
22 65 
18 52 


Annealed 

Annealed and aged 
Cold-drawn 
Cold-drawn and aged.. 


122,200 
101,400 
145,500 























machined with %/4 in. X10 threads. The reduced 
section was finished by grinding and then polished 
in the axial direction, using No. 400 Alumina polishing 
paper. Before testing, each creep specimen, cold- 
drawn or cold-drawn and aged, was heated in a tube 
furnace to the desired temperature and equilibrated 
at temperature for 16 hr. It was then tested at tem- 
perature under constant load. The creep-rupture 
machines used were of a multiple lever type with 
motorized jacks to prevent shock loading. They were 
recently calibrated at room temperature and the 
errors for the stresses used were less than 0.5 percent. 
The temperature of each of the specimens was con- 
trolled to an accuracy of +2 °F. Strain-time data 
were obtained from electric contact follow-up type 
extensometers attached to shoulders of the specimens. 
The sensitivity of extension measurement was ap- 
proximately 0.00005 in. 


3. Metallography 


Most of the specimens tested to fracture in creep 
were subjected to a metallographic examination. In 
general, it was observed that transcrystalline cracks 
having a round or elliptical shape were predominant 
at 700 °F. A mixture of transcrystalline and inter- 
crystalline cracks was evident at 900 °F. The inter- 
crystalline cracks were more populous at the slower 
creep rates. At 1200 °F, the cracks were almost 
exclusively intercrystalline. 
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Four specimens, two each of the annealed and aged 
and the cold-drawn and aged metal, were selected for 
a comprehensive metallographic examination after 
being tested to fracture in creep. The microstruc- 
tures are shown in figures 2 to 7. 
the alloy in the annealed and annealed and aged condi- 
tions were previously shown [1]. 

The basis for selection was (1) these specimens 
illustrate the effect of cold-drawing on the aged micro- 


Microstructures of 


structure; (2) they represent the microstructures asso- 
ciated with the extremes of temperatures used (700 
and 1200 °F); (3) they show the effects of slow creep 
rates on changes in n.icrostructure; and (4) they repre- 
sent microstructures that differed significantly from 
those previously shown for fast creep rates [1], 2]. 
Cracks, formed at 700 °F, are illustrated by the photo- 
micrographs in figure 2A and B. Cracks formed at 
1200 °F are shown in figure 2C and D. More cracks 


FIGURE 2. Relation of test conditions to number and distribution of cracks in fractured 
creep specimens. 


Longitudinal sections, at axis, near fracture, unetched. 
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FiGURE 3. Relation of test conditions to structures of Duranickel after fracture. 
Longitudinal sections, at axis, near fracture. Etched in a solution containing 25 ml each of HNOs, glacial 


acetic acid, HC] and water and 3 g of NaOH: X 667. 


Arrows indicate intercrystalline cracks at “a”, “b”’, and “ce”. 
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were observed in the specimens that had been cold- 
worked prior to aging. If it is assumed that a greater 
number of precipitates were initially present in the 
cold-drawn and aged metal than in the annealed and 
aged specimens, these additional precipitates could 
act as sources for the formation of the cracks as illus- 
trated by the models shown in a previous paper [], figs. 
9 and 10]. The effect of cold-drawing prior to aging 
on the nature of the substructures within the grains 
of the fractured specimens is shown in figure 3. _ Inter- 
secting strain markings were evident in the specimens 
tested at 700 °F (fig. 3A and B): however, a greater 
number of curved, wavy markings was observed in the 
cold-drawn and aged metal (fig. 3B), supporting the 
hypothesis that more precipitates occurred in the cold- 


drawn metal. Some of the markings did not traverse 
the grain from one boundary to the other. Precipi- 
tates, subgrain boundaries, and slip bands acted as 
barriers to the continuous motion of dislocations across 
the grains as will be illustrated later (fig. 5A) by elec- 
tron microscopy. A _ fine structure was observed 
within the grains of specimens tested at 1200 °F (fig. 
3C and D). As will be shown later (figs. 6 and 7) some 
of this substructure consists of overaged precipitates. 
Intercrystalline cracks are larger for the cold-drawn 
and aged specimen (fig. 3Da) than for the annealed and 
aged specimen (fig. 3Cb). Smaller cracks within the 
vicinity of the precipitates at the grain boundaries 
(fig. 3Cc) indicate that voids are nucleated first in 
these regions, then coalesce to form holes and cracks. 


198 





Baye 
Ses ¥ 4.0% of 
Sh Pay Pond, ah t & 
Ht ail the 
eee TW AEN CORA: 
PERCE ea 
REA a Nout: aA 
Oe Rtgs 


LALA 


“ 


“~*> 
a, r 
Ps 
~ 
. . 4 
See ‘ bo 4 
.Y 7 oy 
a a 
rs er 
vs 


FIGURE 4. Electron micrographs showing structure of an annealed and aged alloy 
specimen after fracture in 1367.2 hr at 700 °F. 


Longitudinal sections, near fracture; etched in a solution containing 25 ml each of HNOs, glacial acetic acid, 
HCI and water and 3 g of NaOH; X 7,500. Arrows indicate Ni3Al precipitates at the triple points ‘‘a’’ and “‘b”. 


“ce” indicates substructure near grain boundary. 


Electron microscopy replica techniques were used to 
ascertain more fully the size, number and distribution 
of the precipitates, subgrains, and strain markings in 
these specimens. The results are shown in figures 
4 to 7. Microstructures of the annealed and aged 
specimen tested at 700 °F (fig. 4) indicated the presence 
of annealing twins (fig. 4A), well-defined subgrains, 
- precipitates within some of the grains and an absence 
of them in other grains. It is probable that most were 


too small to be resolved at this magnification. 
tates are also evident at the grain boundaries and 
appear to be largest at the triple points (fig. 4Aa and 


Precipi- 


Bb). Due to differences in structure of adjacent 
parent grains, the shape of the subgrains near the 
boundary is different from that in the interior of the 
grains (fig. 4Bc). When the microstructures shown in 
figure 4 are compared with those of the cold-drawn and 
aged metal, shown in figure 5, it appears that testing 
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FIGURE 5. Electron micrographs showing structure of a cold-drawn and aged 
alloy specimen after fracture in 1068 hr at 700 °F. 


Logitudinal sections, near fracture; etched in a solution containing 25 ml each of HNOs, glacial acetic acid, 
HC1 and water and 3 g of NaOH: X 6,000. Arrows indicate Ni,Al precipitates. 
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FIGURE 6. Electron micrographs showing structure of an annealed and aged alloy 
specimen after fracture in 2059.7 hr at 1200 °F. 


Longitudinal sections, near fracture; etched in a solution containing 25 ml each of HNOs, glacial acetic acid, 
HCl ard water and 3 g of NaOH; X 6,000. “a” indicates peniebanotor region near grain boundary. 
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FIGURE 7. 


Electron micrographs showing structure of a cold-drawn and aged 


Duranickel specimen after fracture in 2452.7 hr at 1200 °F. 
Logitudinal sections, near fracture; etched in a solution containing 25 ml each of HNOs, glacial acetic acid, 


HCI and water and 3 g of NaOH; X 6,000. “ce” 


the aged metal at 700 °F increases the width of the 
subgrain boundaries (fig. 5A and C), the size of some of 
the precipitates (fig. 5C and D, arrows) and the num- 
ber of slip bands (fig. 5B). The discontinuous nature 
of some of the strain markings is illustrated in figure 
5A. The differences in microstructures had a marked 
effect on the creep behavior at 700 °F as indicated by 
the data in figures 8 to 11. 

In general, it was found that the size of the precipi- 
tates increased with increase in temperature and with 
increase in time. The size of the precipitates, ob- 
served in specimens tested at 1200 °F, was less for the 
cold-drawn and aged (fig. 7) than for the annealed and 
aged specimens (fig. 6). Some areas adjacent to the 
grain boundaries (fig. 6Aa) or twin boundaries (fig. 
7Cc) are comparatively free of precipitates. However, 
precipitation within the grains appears to have taken 
place on definite preferred planes and the size of 
these precipitates was less affected by cold-drawing 
than the size of those at the grain boundary. These 
differences in microstructures had little effect on the 
creep behavior at 1200 °F as indicated by figures 9 
and 10. 


4. Test Data and Discussion 


Davies, Richards, and Wilshire [3] studied the effect 
of prior cold-work on creep and fracture of a nickel-0.1 
atomic percent gold alloy. They concluded that, 
below the recrystallization temperature, creep rate 
was reduced and rupture life was increased when the 
prestrain exceeded the initial strain. Their results 


indicates precipitate-free region near twin boundary. 


were interpreted on the basis of varying contributions 
from the grain interiors and the grain boundary regions. 
Dennison and Wilshire [4] carried out tests on three 
grades of nickel and observed that decrease in purity 
was accompanied by an increase in creep resistance 
and rupture life. They concluded that impurities 
impeded grain boundary migration, caused a decrease 
in ductility and an increase in the occurrence of crack- 
ing. The effects of cold-drawing on creep of some 
nickel-copper alloys were discussed by Jenkins and 
Willard [5]. 

Strain-time curves for some of the creep specimens 
used in the present investigation are shown in figure 8. 
Cold-drawing has reduced the extent (time and strain) 
of both the first and second stages of creep below the 
values for the annealed specimens, previously reported 
[2]. For the aged material the first stage was prac- 
tically nonexistent. The shape of the curves are 
both temperature and stress dependent. Time lag 
before the initiation of plastic deformation at 1200 °F 
is believed to be due in part to the blocking of the 
motion of dislocations by the precipitates that formed 
at 1200 °F before loading the specimens. These 
precipitates have been shown both by x-ray and elec- 
tron diffraction to be predominantly particles of Ni3Al. 
The properties of this compound have been studied 
extensively by Guard and Westbrook [6]. 

The relation between stress and second stage creep 
rate? is shown in figure 9. A linear relationship be- 


’The second stage creep rates were obtained from strain-time curves when the data were 
plotted on a linear scale (not shown). 
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FIGURE 8. Strain-time relations of cold-drawn and cold-drawn and aged alloy specimens tested at different 
2 temperatures. 


tween these variables was not obtained when the data 
were plotted on linear, log-log, or semilog coordinates. 
The sigmoidal form of some of the curves indicates the 
existence of a creep rate for which strain hardening is 
a maximum. Cold-drawing the metal increased the 
strength at 700 and 900 °F with the greatest increase 
occurring at the faster rates (fig. 9A). Furthermore, 
cold-drawing prior to aging improved the creep strength 
(fig. 9B) in comparison to the annealed and aged metal. 
No significant difference in strength is observed at 
1200 °F where recrystallization and overaging occurred 
during the tests. 

The effect of stress on time to fracture is shown in 
figure 10. Obviously, no simple equation could be 
used to describe the relationship between these vari- 
ables or between the family of curves for specimens 
having the same or different heat-treatments. It is 
observed, at 700 and 900 °F, that strengthening due 
to cold-drawing or aging is less effective at the longer 
test times than at the shorter times. Recovery of 
the strained material is believed to be the major cause 
of this phenomenon. Although the stress-creep rate 


and stress-fracture time curves appear to be similar, no 
consistent relation existed between creep rates and 
fracture times for either the cold-drawn or cold-drawn 
and aged specimens. 

The interrelationship between second stage creep 
rate and ductility is shown in figure 11. In general, 
elongation and reduction of area values increased with 
increase in creep rate at 700 °F and decreased with 
increase in creep rate at 1200 °F. A large number of 
mechanisms, such as strain hardening, recovery, 
aging, transcrystalline and intercrystalline cracking, 
accompanied creep at 900 °F and caused the ductility 
values to vary in the same manner as was observed 
for annealed and annealed and aged Duranickel [2]. 

Rockwell hardness tests were made at room tem- 
perature on specimens after fracture to ascertain the 
effect of creep on hardening. Hardness indentations 
were made on cross sections of both the shoulder and 
on the deformed (reduced) section near fracture of 
the cold-drawn and cold-drawn and aged specimens in 
order to ascertain the effect of plastic deformation on 
aging characteristics. Small  load-free control 
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FicuRE 9. Effect of cold-drawing and aging on the stress-second stage creep rate relations of the alloy. 
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FIGURE 10. Effect of cold-drawing and aging on the stress-fracture time relations of the alloy. 
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Influence of aging on the ducility-second stage creep rate relations of cold-drawn alloy specimens 


tested to fracture at different temperatures. 


samples, heat-treated for times corresponding to the 
various fracture times in the creep tests, were tested 
for hardness and used for comparison purposes. 
As indicated by the hardness data for the cold-drawn 
specimens (fig. 12A, C, and E) precipitation-hardening 
was more predominant at 900 than at 700 °F. This 
was evident for the specimens exposed te tempera- 
ture only (fig. 12A), for unstrained portions of speci- 
mens (shoulders) exposed to temperature and stress 
(fig. 12C), and for the strained portions of the speci- 
mens (reduced sections) tested to fracture in creep 
(fig. 12E). Extensive recovery, recrystallization, 
and overaging occurred in the specimens tested at 
1200 °F as shown by the decrease in hardness with 
increase in test time. This behavior was also evident 
for the cold-drawn and aged specimens (fig. 12B, 
D, and F) tested at 1200 °F; however, test time had 
only a slight effect on hardness values of the specimens 
tested to fracture at 700 and 900 °F. These data 
tend to confirm the general trends associated with the 
hardness values obtained for the load-free specimens 
and shown in figure 1. 


5. Summary 


Creep tests were made at temperatures of 700, 900, 
and 1200 °F (644, 755, and 921 °K) on a nickel-4.2 per- 
cent aluminum alloy initially as cold-drawn and as 
cold-drawn and age-hardened. 

Microstructures and test data obtained in the present 
investigation are compared to results shown previ- 
ously for the annealed, and annealed and aged alloy. 

The hardness of the alloy was increased by cold- 
drawing or aging. Additional hardening was induced 
by aging subsequent to the cold-drawing operation. 

Due to the existence and formation of numerous 
barriers to the motion of dislocations through the 
crystal lattice, cold-drawing and aging increased the 
creep strength and decreased the ductility of the alloy 
at 700 and 900 °F. No significant effect was apparent 
at 1200 °F, which is slightly above the recrystallization 
temperature. 

Elongation and reduction of area values increased 
with increase in creep rate at 700 °F and with decrease 
in creep rate at 1200 °F. 
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FIGURE 12. Effect of test temperature and time on hardness, at room temperature, of cold-drawn and 


cold-drawn and aged alloy specimens. 


A, B, load-free control specimens as heat-treated; C, D, E, F, creep specimens after fracture. 


Aging the cold-drawn alloy made it more ductile at 
1200 °F and less at 700 °F. 

Post-test hardness values decreased with increase 
in test time at 1200 °F, due to recrystallization. With 
few exceptions, the hardness after testing at 900 °F was 
equal to or greater than that after testing at 700 °F. 

The tendency toward intercrystalline cracking in- 
creased with increase in test temperature. The num- 
ber of cracks was increased by cold-drawing prior to 
age-hardening. 

Microstructures, associated with high strength and 
ductility at the lower temperatures indicated the exist- 
ence of small precipitated particles, small subgrains 
and rather straight intersecting and nonintersecting 
strain markings. However, at 1200 °F, few markings 
were evident and coalesced precipitates appeared 
predominantly in the grain boundaries. The shape 
and number of the precipitates appeared to depend 
on where they are located in the crystal lattice. 


(Paper 70C3-231) 


The authors gratefully acknowledge the cooperation 
of James Heywood and David Ballard in carrying out 
the electron microscopic examination of the specimens. 
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Magnetic Transformation and the Influence of Plastic Strain 
on the Shear Modulus of Fe-Cr-Ni alloys 
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The Fe-18Cr-9Ni and Fe-19Cr-9Ni alloys exhibit a decrease in shear modulus between 76 and 


20 °K. 


occurs at about 40 °K in these alloys. 


Susceptibility measurements confirm that a paramagnetic to antiferromagnetic transition 
The shear modulus for the Fe-25Cr-21Ni alloy did not decrease 


between 76 and 20 °K which was consistent with the fact that the alloy remained paramagnetic to 
4°K. The effect of cold-work on the shear modulus above and below the Néel temperature is discussed. 


Key Words: Antiferromagnetic, cold-work, Fe-Ni-Cr alloys, Néel temperature, shear modulus, 


stainless steel. 


1. Introduction 


It has been shown that the shear and Young’s 
modulus of a Fe-19Cr-9Ni (weight percent) steel 
exhibits a ‘decrease between 76 and 20 °K [1].! Re- 
cently, investigators [2, 3, 4] have reported a para- 
magnetic to antiferromagnetic transition in certain 
18Cr-9Ni steels at about 40 °K. Since it is expected 
that such a magnetic transition would affect the elastic 
moduli, a study was made to better understand the 
effect of the transition on the shear modulus of Fe- 
Ni-Cr alloys. These alloys were Fe-18Cr-9Ni, Fe- 
19Cr-9Ni, and ~~ Fe-25Cr-21Ni_ and _— correspond, 
respectively, to A.I.S.I. commercial stainless steel 
grades 303, 302, and 310. In addition, the effect of 
prior plastic strain at ambient temperatures on the 
change of shear modulus in the temperature range 
295 to 20 °K was documented. 

Some data have been reported concerning the com- 
bined effects of a magnetic transition and _ plastic 
strain on the decrease of elastic moduli on cooling 
below the transition temperature. Various workers 
have noted a change in the Young’s and shear modulus 
of materials that have been cooled at zero magnetic 
field through the ferromagnetic Curie temperature 
[5-10]. Others have noticed the effect of a paramag- 
netic to antiferromagnetic transition on the elastic 
properties of Cr [1]l, 14, 15] and Cu-Mn alloys [12]. 
Recently, Roberson and Lipsitt have reported that 
the Young’s modulus and the shear modulus behave 
differently when Cr is cooled through and below the 
Néel temperature [16]. They found that the AE 
(the change in modulus on cooling through the transi- 
tion temperature) for the shear modulus was less than 


*Cryogenics Division, National Bureau of Standards, Boulder, Colo, 
' Figures in brackets indicate the literature references at the end of this paper. 


217-689 O-66—4 


for the Young’s modulus. Several investigators have 
noted that the AE effect is less for cold-worked ma- 
terials than annealed materials in Fe-Ni-Mo [6], Cr 
[16, 18, 19], and Ni [17]. 

Many have noted a change in elastic modulus after 
plastic strain at ambient temperatures. Although 
some anomalies exist, results on copper, [22, 23, 24, 
27-30] aluminum, [20, 21, 25, 26] and silver [27] in- 
dicate that the elastic moduli are lowered by as much 
as 15 percent when the material is previously plas- 
tically strained. 

Susceptibility measurements described in this paper 
show that the Fe-18Cr-9Ni and Fe-19Cr-9Ni alloys 
become antiferromagnetic at 43 and 39 °K, respec- 
tively, and that the Fe-25Cr-21Ni alloy remains para- 
magnetic down to 4°K. The magnitude of decrease 
in modulus between 76 and 20 °K for the 18Cr-9Ni and 
19Cr-9Ni alloys was found to be inversely related to 
the amount of prior cold-work. 


2. Procedure and Results of Susceptibility 
Measurements 


The susceptibility samples were made from the same 
bar stock that was used for the modulus specimens. 
The composition, hardness, condition and grain size 
of the alloys are listed in table I. The susceptibility 
was measured by the ballistic method. Values of 
susceptibility are plotted in arbitrary units; the units 
of flux change can be converted to susceptibility values 
through calibration of the apparatus A continuous 
temperature change of the sample was achieved by 
allowing the sample to warm from a liquid helium bath. 
The temperature gradient was less than | °K in the 
sample upon warming. Values of susceptibility were 
reproducible with +1 percent: the temperature was 
reproducible within + 1 °K. 
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Table I. Identification of Alloys Tested 


Chemical Composition (weight %) 


Rockwell 
Hardness 





Alloy y S Ni Cr Mo 


Mn Si 


Cu (B scale) Condition 





19Cr-9Ni 
(Commercial 
AISI 302) 


0.02 0.01 8.7 18.6 


0.6 0.6 R, 90 


Annealed 
(~0. 5% martensite) 


Drawn 14% at 
room temperature 
(~13% martensite) 


B 





18Cr-9Ni 0.10 
(Commercial 
AISI 303) 


0.03 0.29 8.7 


17.6 0.4 1.2 0.6 0.4 


Annealed 


Drawn 5% at 
room temperature 


Drawn 13% at 
room temperature 





25Cr-21Ni 
(Commercial 
AISI 310) 


0,08 0:02 0.02 20:8 24:8 0.12 31.7 0.7 0.1 


Annealed 


Drawn 11% at 
room temperature 





Results of the susceptibility measurements for 19Cr- 
ONi, 18Cr-9Ni and 25Cr-21Ni alloys are shown in 
figure 1. The annealed 19Cr-9Ni and annealed 18Cr- 
QNi were measured with a field of about 1 kG; an- 
nealed and 11 percent cold-drawn 25Cr-21Ni were 
measured with a field of about 0.4kG. The 19Cr-9Ni 
specimen, when plastically strained 14 percent, ex- 
hibited a partial phase transformation to body-centered 
cubic martensite. The body-centered cubic phase is 
ferromagnetic [31]. Therefore, the susceptibility 
measurements to determine the antiferromagnetic 
transition were thwarted by the stronger ferromag- 
netism from the martensite phase. Although the 
18Cr-9Ni alloy has similar Cr and Ni content, the 
amount of other stabilizing elements (see table I) is 
greater, which explains the lack of detectable amounts 
of martensite formed by plastic deformation. Since 
the susceptibility versus temperature measurements 
for the 18Cr-9Ni cold-drawn specimens were ap- 
preciably smeared out, it was difficult to detect a dis- 
tinct inflection temperature. 


3. Procedure and Results of Shear Modulus 
Measurements 


The method of determining the shear modulus has 
been previously described [32]. Briefly, the method 
consisted of applying small incremental torque loads 
to carefully machined specimens and measuring the 





19Cr-9Ni , Annealed 


o @ 


FLUX CHANGE, arbitrary units 
S 


25Cr-2INi , Annealed 











40 
TEMPERATURE ,°K 
FicurE 1. The flux change, proportional to the magnetic suscepti- 


bility, as a function of temperature for 19Cr-9Ni steel, 18Cr-9Ni 
steel, and 25Cr-21Ni steel. 
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angle of displacement using an optical lever. The 
maximum shear stress applied to all the samples was 
2900 psi. Measurements were made at three tempera- 
tures: 295, 76, and 20°K. Results are shown in figures 
2, 3, and 4. 
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FIGURE 2. The shear moduli of 19Cr-9Ni, 18Cr-9Ni, and 25Cr-21Ni 


steels as a function of temperature. 
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FIGURE 3. The sheav moduli of 18Cr-9Ni and 25Cr-21Ni steels at 


constant temperature as a function of prior plastic strain. 
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FIGURE 4. The shear modulus of 19Cr-9Ni steel at constant tempera- 
ture as a function of prior plastic strain. 


For any one sample and at any one temperature, the 
deviation from the average value of modulus is 0.5 
percent. The maximum absolute error, i.e., the error 
that occurs either on repeating measurements after 
removal and remounting of the sample, or on measur- 
ing another specimen of the same alloy, is + 2 percent. 


4. Discussion 


In figure 1 it is seen that the 19Cr-9Ni and 18Cr-9Ni 
alloys each show a peak in the susceptibility-tempera- 
ture curve which is characteristic of a material trans- 
forming from a paramagnetic to an antiferromagnetic 
state. These results are consistent with the modulus 
curves of figure 2; the 19Cr-9Ni and 18Cr-9Ni alloys 
exhibit anomalous drops in modulus whereas the 25Cr- 
21Ni alloy behaves normally down to 20 °K. The de- 
crease in modulus between 76 and 20 °K is due to an 
additional strain component resulting from antiferro- 
magnetic ordering below the Néel temperature [12, 13]. 

From the 18Cr-9Ni and 19Cr-9Ni alloy shear modulus 
curves of figure 2, it is seen that (Gze— Goe)/Gre is 


- inversely related to the amount of prior plastic strain: 


the AG/G for the 18Cr-9Ni alloy is respectively 4.9, 
3.6, and 2.6 percent for annealed, 5 percent cold-drawn, 
and 13 percent cold-drawn conditions; AG/G for the 
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19Cr-9Ni alloy decreased from 3.6 to 1.4 percent for 
the annealed to cold-drawn conditions. This decrease 
indicates that there is less magnetostrictive effect 
at 20 °K in a cold-worked sample due to internal 
stresses from defects which influence the redistribu- 
tion of the domain vectors [35, 6]. 

The effect of prior plastic strain above the Néel 
temperature on the shear modulus is illustrated in 
figures 2to 4. A decrease of the modulus values, after 
straining, has been verified. This is expected since 
many of the dislocations introduced during plastic 
straining move under very low stresses. In addition, 
vacancies created by cold-work lower the material 
density which lowers the modulus. Jongenburger 
[33] has predicted that materials with low stacking 
fault energy would exhibit greater sensitivity to cold- 
work than high stacking fault energy materials, since 
the stress needed to move extended dislocations should 
be less. The 18Cr-9Ni alloy is expected to have a 
lower stacking fault energy than the 25Cr-21Ni alloy 
[34]. At ambient temperatures the strain sensitivity of 
the shear modulus for the 18Cr-9Ni alloy is greater: 
however, at 76 °K their sensitivities are approximately 
equal. The reason for this temperature trend is not 
clear. + 

At 20 °K, notice that the shear modulus of the 18Cr- 
QNi alloy is essentially independent of the amount of 
prior plastic strain. Apparently the AE effect and the 
defect contribution, which serves to lower the overall 
modulus, just offset each other at this temperature. 
The much greater reduction of the shear modulus 
of the 19Cr-9Ni alloy as a function of strain is un- 
doubtedly influenced by the partial transformation to 
the martensite (body centered cubic) phase during 
plastic deformation at room temperature. 


5. Summary 


It has been found that the elastic properties of low 


stacking fault energy austenitic stainless steels 
undergo a transition at about 40 °K. This transition 
preduces a decrease in the shear modulus values and 
is caused by a paramagnetic to antiferromagnetic 
transformation. The decrease of the modulus values 
between 76 and 20 °K is inversely related to the amount 
of prior cold-wogg at ambient temperature. The 
shear modulus at temperatures above the Néel tem- 
perature is influenced by previous cold-work, decreas- 
ing up to 10 percent for 10 to 15 percent plastic strain. 


The authors thank J. J. Gniewek for informal dis- 
cussions regarding the work in this paper. 
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Thermal relaxation and Brillouin scattering in liquids, R. D. 
Mountain, J. Res. NBS 70A3 (Phys. and Chem.), No. 3, 207-220 
(May-June 1966). 

Weak coupling of internal degrees of freedom of molecules to 
the translational degrees of freedom of a fluid results in additional 
modes of motion for density fluctuations. These new modes 
affect the spectral distribution of light scattered by density fluctua- 
tions so that the Landau-Placzek ratio is not satisfied. The case 
of thermal relaxation with a single relaxation time is worked out 
in detail. Formulas for the spectral distribution of the scattered 
light, for the ratio of the intensities of the central (Rayleigh) to 
the Brillouin components and for the phonon velocity are derived 
and applied to carbon disulfide and carbon tetrachloride. The 
data for carbon tetrachloride are shown to be inconsistent with 
the single relaxation time model for thermal relaxation. 


Calibration of germanium resistors at low temperatures 
(2-20 °Kelvin), G. Cataland and H. H. Plumb, J. Res. NBS 70A3 
(Phys, and Chem.), No. 3, 243-252 (May-June 1966). 

Two germanium resistors, C and D, have been referenced to values of 
temperature determined by the NBS acoustical thermometer. The 
resistance-temperature calibration data for these resistors have been 


fitted to a function of the form logo R= A,(logioT’)" and the 
n=0 

results are presented. The resistors C and D that are used as 
secondary standards maintain the scale, NBS Provisional Scale 2 
to 20 (1965), to which public calibrations of germanium thermometers 
are referenced. 

The calibration apparatus and measurement techniques that are 
employed in calibrating submitted resistors are described in detail. 
Additionally, data of three typical calibrations and their treatment 
by polynomial fitting are included to demonstrate the characteristics 
of some commercially available germanium resistors. 


Mechanical relaxation in polyethylene crystallized with 
various degrees of lamellar orientation, J. M. Crissman and 
E. Passaglia, J. Res. NBS 70A3 (Phys. and Chem.), No. 3, 225- 
232 (May—June 1966). 

The mechanical relaxation behavior of a set of well-characterized 
samples of polyethylene crystallized with different degrees of 
lamellar orientation is reported. The various samples ranged in 
morphology from unoriented isotropic samples to ones which showed 
a high degree of orientation of the b-axis along the sample growth 
direction. The mechanical measurements were made using a 
torsion pendulum apparatus of standard design, the direction of shear 
being normal to the b-axis for the oriented samples. The tempera- 
ture range covered was from 100 to 400 °K. No definite effects 
attributable to orientation were observed for either the y or B 
relaxation process, whereas for the a@ relaxation results for G’ 
indicate that a slight decrease in peak height resulted from the 
presence of lamellar orientation, particularly on the high tempera- 
ture side of this peak. Data for the real and imaginary parts of 
the complex shear compliance are also discussed. 


On a sequence of points of interest for numerical quadra- 
ture, S. Haber, J. Res. NBS 70B2 (Math. and Math. Phys.), No. 2, 
127-136 (Apr.—June 1966). 

The distribution of 2 sequences of points, originally discussed by 
van der Corput and by K. F. Roth, is studied in detail. The results 
obtained disprove a conjecture of J. H. Halton and suggest a con- 
jecture on the improvement of a theorem of Roth. 


Treatment of outliers in samples of size three, F. J. Anscombe 
and B. A. Barron, J. Res. NBS 70B2 (Math. and Math. Phys.), 
Vo. 2, 141-147 (Apr.—June 1966). 
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Let X be a space of functions, say X C C(K), K locally compact 
Hausdorff, let /eX* be an integral on X and let M* C X* be a given 
subspace of “simple” functionals, then it is desired to obtain an 
leM* for given n, leM* C M*; M* being a suitable n dimensional 
subspace determined so that / if /— / annihilates a given finite dimen- 
sional subspace X,CX. In this general context, the abstract 
analyses of numerical intergration is developed and certain specific 
applications are made. 


Radio meteorology, B. R. Bean and E. J. Dutton, NBS Mono. 92 
(Mar. 1, 1966). $2.75. 

This volume brings together the work done in radio meteorology 
over the past ten years at the National Bureau of Standards’ Central 
Radio Propagation Laboratory (CRPL). This decade has seen the 
development, on an international scale, of great emphasis upon the 
effects of the lower atmosphere on the propagation of radio waves. 
The CRPL group has concentrated upon the refraction of radio 
waves as well as the refractive index structure of the lower atmos- 
shere on both synoptic and climatic scales. These are the areas 
of radio meteorology that are treated in this volume, with additional 
chapters on radio-meteorological parameters and the absorption of 
radio waves by the various constituents of the lower atmosphere. 


Copper wire tables, NBS Handb. 100 (Feb. 21, 1966). 
Supersedes Circ. 31. 

This Handbook is a revision of the Copper Wire Tables previously 
published as NBS Circular 31. It reflects changes in the nominal 
diameters of gages 45 and smaller and extends the tables to 56 gage. 
The changed diameters and extended range were established in 1961 
by the Committee on Wires for Electrical Conductors of the Ameri- 
can Society for Testing and Materials and were published as ASTM 
Standard B258-61. They have also been approved as American 
Standard C7.36-1961 by the American Standards Association. To 
reduce internal inconsistencies, tables 5 through 14 were completely 
recomputed by the ASTM Committee on Wires for Electrical 
Conductors. 

The first edition of Circular 31 was published in 1912 at the request 
of the American Institute of Electrical Engineers. Subsequent 
editions appeared in 1914 and 1956. 


50 cents. 


OMNITAB. Acomputer program for statistical and numeri- 
eal analysis, J. Hilsenrath, G. G. Ziegler, C. G. Messina, P. J. 
Walsh, and R. J. Herbold, NBS Handb. 101 (Mar. 4, 1966). $3.00. 
This user’s manual gives the characteristics and application of 
a general-purpose interpretive program for the generation of 
physical and mathematical tables on the digital computer using 
FORTRAN. The vocabulary and sentence structure permit (with 
only a few simple rules) the translation of a mathematical state- 
ment into text understandable both to the user (not necessarily a 
programmer) and to the machine. The basic vocabulary of slightly 
over 250 English words is given. These cover the arithmetic and 
mathematical operations such as ADD, DIVIDE, MULT, SUB- 
TRACT, LOGE. EXP: the elementary and special mathematical 
functions: SIN, TAN, COSH, TANH, orthogonal polynomials, 


Bessel functions, elliptic integrals, etc; and various manipulative 


_operations on columns of numbers and on entire arrays. 


Detailed instructions are given also for applying the program 
to: matrix calculations, statistical analysis of experimental data, 
least-squares curve fitting, linear and nonlinear interpolation, nu- 
merical differentiation and integration, solutions of systems of 
equations, and location of zeros and maxima and minima of algebraic 
and transcendental functions. 

Numerous examples are given of applications of the program to 
scientific calculations covering a wide spectrum of problems en- 
countered in a research and development organization. Aside 
from the comprehensive repertoire of mathematical functions and 
operations, the important features which are stressed are the flexible 
format-free input and output provisions and curve plotting facility. 





Procedures for precise determination of thermal radiation 
properties, November 1963 to October 1964, J. C. Richmond, 
S. T. Dunn, D. P. DeWitt, and W. D. Hayes, Jr., NBS Tech. Note 
267 (Dec. 17, 1965). 45 cents. 

A laser-source integrating sphere reflectometer was designed and 
built to measure the reflectance of specimens at high temperatures. 
The sphere was calibrated for linearity of response at 0.6324 by 
means of a shallow cylindrical cavity with a variable depth-to-radius 
ratio, having a lining of known reflectance. Preliminary tests 
showed that the flux emitted by a hot specimen at temperatures up 
to 2500 °K will not invalidate the reflectance measurement. An 
ellipsoidal mirror reflectometer was calibrated for all known errors 
in the 1 to 7.54 range. A preliminary analysis indicates that the 
errors in measurement of absolute reflectance with this equipment 
should not exceed 2 percent. A review of the literature on relation 
of thermal radiation properties to other properties of materials is 
presented, together with a summary of the work done in an effort to 
compute the spectral emittance of rhodium. 


Earth’s field static calibrator for accelerometers, P. S. 
Lederer and J. S. Hilten, NBS Tech. Note 269 (Feb. 1, 1966). 20 
cents. 

This paper describes a simple, relatively inexpensive system for 
the precise static calibration of accelerometers in the earth’s gravi- 
tational field with an estimated limit of error of +0.0004 g. The 
system consists of a precision machinists dividing head on a surface 
plate and a precision level. Experimental results obtained during 
the calibtation of three types of accelerometers are shown to indi- 
cate the capabilities and usefulness of the system. 


Radiochemical analysis: Activation analysis, instrumenta- 
tion, radiation techniques, and radioisotope techniques, 
July 1964 to June 1965, ed. J. R. DeVoe, NBS Tech. Note 276 
(Jan. 7, 1966), $1.00. 

This is the second summary of progress of the Radiochemical 
Analysis Section of the Analytical Chemistry Division at the Na- 
tional Bureau of Standards. 

Arrangements for the new facilities at the reactor building and 
Linac are in their final stages of completion. Additional work is 
described on the electronics for the handling of data from pulse 
height analyzers as well as the use of digital computers for processing 
these data. 

A Mossbauer spectrometer which uses optical interferometry in 
the Doppler drive, and the mechanical design for the electromechan- 
ical drive described in last years report are presented. A descrip- 
tion of chemical shift of iron compounds is given. The use of a 
correction factor for certain source detector geometries are described 
for Méssbauer spectroscopy. A new source of Sn!!9” in palladium 
which produces a very narrow line width is also described. 

For the first time a trace element has been certified in a Standard 
Reference Material with the assistance of the activation analysis 
technique. Flux measurements for the NRL reactor are given. 
A number of special analyses for various NBS projects by activa- 
tion analysis are described (e.g., argon in silicon, chlorine in poly- 
styrene, hafnium in zirconium). A beta-gamma sum coincidence 
gamma-ray spectrometer is described. 

A constant potential electrochemical method of analysis using 
radioisotopic tracers is presented which has potential high sensi- 
tivity. A method of using radioisotopic tracers in heterogeneous 
equilibrium between an aqueous medium of cobalt to be analyzed 
and manganese dioxide impregnated filter paper is described. 
Some results on the precision of the substoichiometric radioisotope 
dilution analysis of silver and iron are presented. 


A peak ac-de voltage comparator for use in a standards 
laboratory, L. A. Marzetta, NBS Tech. Note 280 (Jan. 17, 1966). 
25 cents. 

An increased availability of high-precision a-c and d-c power sources 
has allowed routine electrical measurements to be made in the field 
with a precision previously possible only in the standards lab- 
oratory. Ac-de comparisons are being made on the basis of the 
peak, rms, and average value of the sine wave source. The peak 
ac-de voltage comparator described in this report was developed 
for use in the standards laboratory to help meet the increasing 


demand for the certification of new field instruments. The de- 


scribed instrument is capable of establishing the ac-dc error of 
other peak comparators at frequencies from 50 to 2400 hz with an 
imprecision of less than 15 ppm. 


Keypunch controls for string-punching of statistical data, 
M. J. Brennan, NBS Tech. Note 282 (Feb. 11, 1966). 25 cents. 
This note describes controls for a keypunch machine used to key 
statistical data with “string-punching techniques.” Since the data 
are punched into cards that do not contain fixed field sizes, an 
identification code is punched for each data entry. To avoid a split 
field (part on one card and the remainder on another), the keypunch 
machine is programed to provide a justification zone in which an 
answer identification code or a “delete” code cannot be keyed. 
The cards are automatically serialized for each questionnaire and 
an arrangement is provided to punch an end-of-questionnaire code 
in the last card for the questionnaire. 


An outline for cooperative action on the determination of 
x-ray mass attenuation coefficients in the wavelength range 
from 0.5 to 100 A (25 to 0.12 keV), K. F. J. Heinrich, NBS 
Tech. Note 284 (Feb. 21, 1966). 15 cents. 


. Sources of information on x-ray mass attenuation coefficients and 


the present status of data in this area are reviewed. Many gaps 
are observed in the existing data and some values derived by in- 
terpolation may be of questionable accuracy. Because of the 
importance of accurate mass attenuation coefficients, the extension 
and improvement of these measurements is proposed together with 
a plan for coordination of effort and exchange of information. 


Computation of the permeability and permittivity of a rela- 
tively small ring sample in a toroidal coil, E. G. Johnson, Jr., 
NBS Tech. Note 311 (Jan. 25, 1966). 30 cents. 

The derivation and the FORTRAN machine program are presented 
for a formula which gives the impedance of a partially filled toroid 
coil. This formula shows the relationship of the complex permea- 
bility and permittivity of the ring material, the dimensions of the 
coil and sample, and the frequency of the applied electromagnetic 
field. Limitations of the formula due to proximity of the ring and 
coil, range of permeability and permittivity of the ring, and uni- 
formity of the coil winding are considered. The cost of the program 
in terms of computer time is also considered. The principal value 
of the formula is to accurately determine the permeability and 
permittivity of ferrites. 


A 100 kW 2-25 Me/s distributed amplifier, designed for 
use with 10 kW ionospheric sounders, W. B. Harding, M. W. 
Woodward, and J. C. Carroll, NBS Tech. Note 331 (Feb. 15, 1966). 
20 cents. (Formerly Central Radio Propagation Laboratory of 
the National Bureau of Standards, now Environmental Science 
Services Administration, Boulder, Colo.). 

The principles of distributed amplification have been known for 
many years and have been applied to fast rise time pulse and oscil- 
loscope amplifiers for small-signal applications [Percival, 1935-37: 
Pettit and McWhorter, 1961]. This paper describes the design, 
construction, and testing of an r-f distributed amplifier intended 
for ionospheric sounding and backscatter applications. The 
amplifier has a power gain of 10 dB and a rated output of 100 kW 
peak pulse power, over a range of 2 to 25 Mc/s. A feature of this 
amplifier is the elimination of input and output transformers. 


Restoration of complete dentures inadvertently warped by 
the patient: Report of case, J. B. Woelfel and G. C. Paffenbarger, 
J. Am. Dental Assoc. 71, 866-870 (Oct. 1965). 

A set of Hydro-Cast dentures that had given satisfactory service 
for over 5 years was accidently warped when the patient poured hot 


water over them. The warpage was manifested by a molar-to-molar 
shrinkage of 0.43 mm. on the upper and 0.5 mm. on the lower den- 
ture. The corresponding flange-to-flange shrinkages were 1.35 
mm. on the upper and 0.67 mm. on the lower. These changes 
rendered the upper denture unusable. Both dentures were satis- 
factorily altered by selectively grinding areas of the denture base 
which were in hard contact with the supporting tissues. The ap- 
position of the tissues and the dentures was determined with pat- 
terns produced on the tissue-bearing surfaces of the dentures by 
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a thin coating of a pressure-indicator paste. After four trial adjust- 
ments with the upper denture and three with the lower denture, 
they fitted satisfactorily, were comfortable and functioned well 
as evidenced by the patient’s ability to chew such foods as celery 
and to play the trumpet regularly. The supporting tissues had a 
normal appearance after the adjustments. 


Spectrophotometric determination of bromine and hydro- 
gen bromide, E. C. Creitz. Anal. Chem. 37, 1690-1692 (Dec. 1965). 
A method is described for determination of Br. and HBr in the con- 
centration range zero to one microgram Br per ml. The analysis 
depends on oxidation of o-tolidine by Br. and subsequent colori- 
metric determination of the yellow oxidation product. Total 
bromine was determined by measurement of the optical absorption 
of AgBr suspensions, and HBr obtained by difference. An overall 
reproducibility of about 2 absorbance units (approximately 70 nano- 
grams per ml) was obtained. Some oxidizing or reducing sub- 
stances may produce interferences. 


Comment on obscurities of oscillator noise, L. Fey, W. R. 
Atkinson, and J. Newman, Proc. IEEE Letter 52, 104-105 (Jan. 
1964). 

Because of the need for knowledge of oscillator noise character- 
istics and because other published results conflict with those in a 
recent paper by the present authors, a separate analysis of this 
published data was made in order to ascertain whether any alterna- 
tive conclusions might be reasonable. This analysis showed that 
the data might also be interpreted to imply flicker type noise, as 
reported by the present authors for oscillators which they examined, 
rather than noise of a thermal origin as suggested in the original 
publication, 


Cryogenic fluid, two-phase critical flow studies related to 
reactor systems, R. V. Smith, H. S. Isbin, H. K. Fauske, M. Pet- 
ick, C. H. Robbins, and F. R. Zuludek, Proc. 3d United Nations 
Intern. Conf. Peaceful Uses of Atomic Energy, Geneva, Switzerland, 
Aug. 31-Sept. 9, 1964. 

This paper reviews the applications where critical two-phase flow 
of cryogenic fluids is possible in reactor systems and summarizes 
the analytical and experimental work that has been done with 
respect to this problem. 


Field strength calibration techniques at the NBS, H. E. 
Taggart, IEEE Trans. Electromag. Compatability EC-7, No. 2, 
163-169 (June 1965). 

Instrumentation utilized by the National Bureau of Standards to 
calibrate field strength meters and standards in terms of absolute 
field strength at frequencies from 10kHz to 1 GHz are described. 
Two basic techniques are used: (1) the standard-field method for 
loop antennas at frequencies from 30 kHz to 30 MHz and (2) the 
standard-antenna method for horizontally-polarized dipole antennas 
at frequencies from 30 to 1000 MHz. The accuracy and the limi- 
tations of each technique are discussed. Present loop antenna 
calibration uncertainties vary from 3 to 5 percent depending on 
the frequency. Dipole antenna uncertainties are presently 12 per- 
sent. Various other antenna calibration techniques that have been 
studied are mentioned, pointing out their limitations. The instru- 
mentation used to calibrate the field strength receiver character- 
istics is described. These measurements include the calibration 
of the signal attenuators, calibration of the overall linearity of the 


receiver, and calibration of the receiver as a two-terminal rf 


voltmeter. 

Future plans of the NBS field strength and noise calibration services 
are discussed. Improved uncertainties and extended services 
are planned for antennas from 30 Hz to 1000 MHz. New noise 


calibration services are planned; these include the calibration of 


random noise generators at selected frequencies from 30 to 1000 
MHz, and the calibration of impulse generators in terms of their 
spectral density. 


Field-strength measurements in a multipath field, C. C. 
Watterson, Proc. 4th Natl. IEEE Symp. Radio Frequency Inter- 
ference, San Francisco, Calif., June 28-29, 1962. 

Field-strength measurements made in the VHF—UHF portion of 
the spectrum often have limited accuracy and are ambiguous be- 
cause of the prevalence of multipath propagation. This paper 


describes a method of making field-strength measurements in a 
multipath field from which it is possible to calculate the amplitude, 
polarization, relative time phase, and direction of arrival of each 
of the multipath components. The method is basically an aperture- 
synthesis technique which requires that the amplitude and phase 
of the signal from a moving antenna be measured as the antenna 
is moved successively along three orthogonal straight-line paths. 
Three antenna orientations are used on each of the three paths. 
From the measurements of the signal amplitudes and phases it is 
necessary to solve nine parallel, similar problems, one for each of 
the nine combinations of antenna orientation and measurement path. 
The nine solutions are combined to provide the final solutions for 
each multipath component: The amplitude and relative time phase 
of the horizontal electric field component, the amplitude and time 
phase of the “vertical” electric field component, and the azimuthal 
and elevation angles of arrival. 


Gearing errors as related to alinement techniques of the 
rotary-vane attenuator, W. Larson, 19th Annual IEEE Conference 
and Exhibit on Instrumentation and Measurements, New York, N.Y., 
Oct. 12—15, 1964, Preprint No. 12—5—2—64 (Sept. 1965); IEEE Instr. 
Meas. IM-14, No. 3, 117-123 (Sept. 1965). 

The rapid advancement of microwave techniques in the measure- 
ment of attenaution, power, noise, and impedance has increased 
the need for a precision rf variable attenuator. The rotary-vane 
attenuator appears to be potentially capable for use as a high- 
precision rf attenuator or attenuation standard. An analysis of 
the measured errors of a large number of rotary-vane attenuators 
has shown that the systematic error present is correlated with vane 
misalinement. Also, a study of the gear mechanism used to drive 
the rotating vane has revealed another systematic error for which 
corrections can be made. 

Graphs are used to illustrate the readout scale error in decibels at 
various magnitudes of the drive-gear eccentricity and driven gear 
pitch diameter. An electrical alinement technique utilizing a cor- 
rection equal to the “‘average vane angle error,” which is determined 
from calibration data, is used to aline the rotating vane to an accuracy 
of +0.0012°. Vane alinement by this technique can reduce the 
systematic error of vane misalinement by an order of magnitude 
and decrease the deviation in attenuation at high values of attenua- 
tion as much as two orders of magnitude. 


Optical instrumentation for the biologist, J. R. Meyer-Arendt. 
Appl. Opt. 4, No. 1, 1-9 (Jan. 1965). 
A great variety of optical instruments have been developed for use 


in biology. The microscope and its many special types are part 
of them. In this article, a survey is given of the history. develop- 
ment, current status, and future outlook of general microscopy. 
microspectroscopy, interference microscopy. and 


and 
X-ray microscopy. 


electron 


The NBS cryogenic data center, Proc. X/ Intern. Congress Re- 


frigeration, Munich, Germany, Aug. 28, 1963, pp. 110-114 (1965). 


A brief description is given of the two principal activities of the Cryo- 
genic Data Center, namely: the evaluation and compilation of data 
on low temperature properties of materials; and the acquisition, 
handling and dissemination of technical literature. A program for 
thorough evaluation and the compilation of thermodynamic and 
transport properties of cryogenic fluids and selected solids is men- 
tioned. Four functions of the literature service are outlined and 
typical pictures of the operation included. 


Timing and space navigation with an existing ground based 
system, G. Hefley, R. F. Linfield. and R. H. Doherty, AGARD Navi- 
gation Systems for Aircraft and Space Vehicles, pp. 490-506 (Per- 
gamon Press Inc., New York, N.Y., 1962). 

New time synchronization techniques have been developed with 
accuracies which make feasible the use of an inverse hyperbolic 
system for determining the position of space vehicles. 

The National Bureau of Standards has developed a Loran-C clock 
which provides timing and synchronizing accuracies of one micro- 
second or better and stabilities of 0.1 microseconds or better any- 
where within a Loran-C navigation chain. Synchronization of 
other chains could provide base line lengths equal to an earth’s 
diameter. 
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Having established a common time base, short pulse transmissions 
from space vehicles are measured relative to this base. _Differ- 
ences in. the time of arrival of these pulses at various ground sta- 
tions permit space fixes at great ranges. 

The time synchronization scheme is discussed along with fix ac- 
curacies obtainable as a function of synchronization accuracy and 
system geometry. 


22 ISCC-NBS centroid colors of maximum contrast, K. L. 
Kelly, Color Engr. 3, No. 6, 26—27 (Nov.-Dec. 1965). 

Sets of colors having maximum contrast with one another are often 
needed for color coding. The number of colors required will 
depend upon the application. Twenty-two ISCC-NBS centroid 
colors have been selected and arranged in such an order that max- 
imum contrast is obtained by taking the required number of colors 
from the top of the list. 


Elastic stress-strain relations in perfect elastic fluids, B. 
Bernstein, E. Kearsley, and L. Zapas, Trans. Soc. Rheol. 9, Pt. 1, 
27-39 (1965). 

The theory of perfect elastic fluids deals with the nonequilibrium 
¢hermodynamics of finite viscoelastic deformation. The thermo- 
dynamics of finite elastic strain and of perfect elastic fluids are 
discussed. It is shown that the stress-strain relations for perfect 
elastic fluids may be written in a way which appears formally iden- 
tical to those in elasticity theory, but the form of this dependence 
depends on past history. 


A stable are source of high ultraviolet radiance, C. R. Yokley,. 
Proc. Symp. Solar Radiation Simulation, Los Angeles, Calif., Jan. 
8-20, 1965, pp. 107-110 (1965). 

A high pressure argon arc with an intense continuum has been de- 
veloped as a possible source in ultraviolet radiometry. The wall- 
stabilized, continuously operating arc exhibits radiance temperatures 
of about 5000 °K at 2400 A. Typical operating conditions are cur- 
rents from 50 to 100 amperes and pressures from 5 to 25 atmos- 
pheres. Reproducibility of spectral radiance calibration is about 
3 percent when argon flow, arc current, and chamber pressure are 
regulated within prescribed limits. Short term reproducibility 
over a period of several hours is about | percent in spectral radiance. 
The spectral radiance is approximately proportional to the square 
root of the pressure with an upper stability limit of 375 psia for 50- 
ampere arcs. The are design, stabilizing and operating techniques 
along with typical spectral radiance calibrations from 2200 A to 
7000 A is given. 


Coherent Raman effect in the off-axis Raman resonator, 
D. A. Jennings and H. Takuma, Appl. Phys. Letters 4, No. 11, 
185-186 (June 1964). 

The ruby Q-switched laser has been used to excite coherent Raman 
oscillation in a separate optical resonator. The Raman effect was 
excited by the 20 MW ruby laser through a hole in one of the mirrors 
of the Raman resonator. The Raman resonator was placed at a 
small angle relative to the ruby laser beam axis. Experimental 
evidence is given to show that the Raman output is a build-up of 
oscillations along the axis of the Raman resonator. 


Estimation of the parameters of the Hyper-Poisson distri- 
bution, E. L. Crow and G. E. Bardwell, (Proc. Symp. Classical 
and Contagious Discrete Distribution, Montreal, Canada, Aug. 
15-20, 1963), Chapt. V, pp. 127-140 (Statistical Publ. Soc., Cal- 
cutta, India, 1965). 

The hyper-Poisson family of discrete distributions, a two-parameter 
generalization of the Poisson distributions, was introduced recently 
by the authors (Journal of the American Statistical Association, 
59, March 1964). In this paper six different pairs of estimators 
(including four proposed in the previous paper) for the two param- 
eters are presented and discussed: maximum likelihood, two 
moments, three moments, two moments plus zeroth frequency, 
mean plus first two frequencies, and two frequencies only. Three 
of these pairs are given by simple explicit formulae. Asymptotic 
variances are obtained for the three-moment estimators, and em- 
pirical variances for the three pairs of explicit estimators. Some 
of the estimators are applied to five examples of experimental or 
survey data. 


Generalized variational principles for electromagnetic 
junctions. Vibrations: application to the theory of wave- 
guide, Proc. Symp. Electromagnetic Theory and Antennas, Copen- 
hagen, Denmark, June 1962, ed. E. C. Jordan, pp. 253-260 (Per- 
gamon Press Inc., New York, N.Y., 1963). 

Generalized variational principles are stated for the elements— 
and for the associated electromagnetic fields—of the immittance 
matrices of waveguide junctions. The formulation applies to junc- 
tions containing arbitrary dissipative and anisotropic media, for 
which the underlying differential equations are not self-adjoint; 
and, by the use of bilinear (rather than quadratic) functional forms, 
permits direct treatment of off-diagonal elements of the immit- 
tance matrices. 


Infrasonic observations of the May 16, 1964, voleanic ex- 
plosion on the island of Bali, V. H. Goerke, J. M. Young, and 
R. K. Cook, J. Geophys. Res. 70, No. 24, 6017-6022 (Dec. 15, 1965). 
The volcanic eruption of Mount Agung on the Island of Bali reached 
explosive proportions on May 16, 1963. Perturbations generated 
in the atmosphere traveling at sonic velocities were recorded in 
Boulder, Colo., 14,700 km distant, both via the short great-circle 
path and via the long great-circle path through the antipodes and 
around the world. In Boston, Mass., 16,200 km distant, and in 
Washington, D.C., 16,300 km distant, sound waves were received 
via the short great-circle path. The average transit velocity for 
the short path was 17 m/sec less than for the long path. The effect 
of wind velocities at mesospheric heights on the transit velocities 
is discussed. 


On the use of thermopiles for absolute radiometry in the 
far ultraviolet, R. G. Johnston, and R. P. Madden, Appl. Opt. 4, 
No. 12, 1574-1580 (Dec. 1965). 

Thermopiles used for. absolute radiometry in the far ultraviolet 
spectral region are generally calibrated with visible and near 
infrared radiation. Three possibilities for the failure of the sen- 
sitivy of thermopiles to be independent of wavelength have been 
investigated. It was found that gold-black typical of that used on 
thermopiles for far ultraviolet applications reflectivity scattered 
2 to 4 percent of the incident far ultraviolet radiation. Since the 
same .gold-black scattered a similar fraction of incident visible and 
near infrared radiation, the correction for this effect was small. 
An analysis of surface sensitivity maps indicated that a thermopile 
of 0.04-sec time constant showed no wavelength variation of sen- 
sitivity when either de detection or ac synchronous detection with 
13 cps chopping was used. The most significant wavelength-de- 
pendent phenomenon was found to be the photoelectric effect. 
Ejected electrons carry a significant amount of energy away from 
gold-blacked thermopiles in the spectral region below 1600 A. The 
maximum correction determined for a particular thermopile was 
5.2 percent at A 735 A. 


Solute nuclear magnetic resonances in primary lead alloys, 
L. H. Bennett, R. M. Cotts, and R. J. Snodgrass (Proc. XIIIth 
Colloque and Relaxation in Solids), Book, Nuclear Magnetic Res- 
onance and Relaxation in Solids, ed. L. V. Gerven, pp. 171-180 
(North Holland Publ. Co., Amsterdam, The Netherlands, 1965). 
The nuclear magnetic resonances of In' Bi2, T2%, Hg!, and 
Sb!*! as solutes in fee primary lead solid solutions have been meas- 
ured. The room temperature Knight shifts for these solutes at 
low concentration in Pb are: In' >, k=0.86%; Bi?, k=1.15%; 
TP®, k= 1.93%; Hg, k= 2.41% and Sb'?!, k= 0.66%. In all cases 
the change of Knight shift with temperature, magnetic field, or con- 
centration, if any, is small. Quadrupole effects, if any, are small. 
From the temperature dependence of the linewidth for Tl?%, 7, is 
found to be about 7 w sec, compared with 7.3 yw sec found by use 
of the Korringa relation. Comparison of the solute shifts is made 
both with the solvent shift, and with the shift of the solute in itself, 
assuming a free electron model. 


Brush cathode plasma—a well behaved plasma, K-B. Persson, 
J. Appl. Phys. 36, No. 10, 3086-3094 (Oct. 1965). 

Results of probe, spectroscopic and microwave studies are used 
to describe properties of a novel cold cathode discharge. The use 
of a brush cathode greatly exaggerates some of the usual cold cath- 
ode mechanisms, thus determining the interesting features of 
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the plasma produced. ‘The cathode fall is about an order of magni- 
tude larger than for the corresponding normal cold cathode. The 
brush cathode generates a uniform electron beam in the energy range 
1 to 10 kilovolts and a corresponding negative glow with a longitudinal 
dimension (reaching distance of the beam) one to two orders of 
magnitude larger than for the normal cathode. 
sions of the negative glow and the fact that it is extremely well- 
behaved (no instabilities and no striations) make it ideal for a whole 
series of investigations in plasma physics and spectroscopy. The 
electron density in the helium plasma is in the range 10!° to 10'4 
cm~*, the electron temperature in the range 0.05 to 0.10 eV. The 
negative glow is beam generated (essentially field free) and recom- 
bination dominated making it a practically uniform plasma ideal 
for the study of rate and transport coefficients. Preliminary 
measurements of the particle balance give a recombination rate of 
5 X 10-'° cm® sec~! at 1100 °K and an electron density of 3 X 10" 
cm? in excellent agreement with the collisional-radiative recombi- 
nation theory developed by Bates, Kingston, and McWhirter. 
The recombination light emitted by the plasma is sufficiently strong 
to make spectroscopic methods very useful for measurements of 
the electron density and the electron temperature. Lines in the 
series 2s°S—np*P® are observed up to the quantum level n=30. 


The realization and use of the international practical tem- 
perature scale of 1948, J. P. Evans, ISA 20th Annual Confer- 
ence and Exhibit, Los Angeles, Calif., Oct. 4—7, 1965), ISA Preprint 
14.3—1-65 (1965). 

This paper describes the International Practical Temperature 
Scale of 1948, the methods and accuracy with which it is realized, 
and the current efforts being made to improve the scale. Methods 
of transferring the scale to temperature measuring instruments 
are also discussed. 


Topside sounding as a tool for global ionosphere studies, 
R. W. Knecht, JEEE Trans. Ant. Prop. Programs Digest, pp. 21-23 
(1963). 

With the advent of the successful Alouette topside sounding satel- 
lite, a powerful new technigue has become available for exploration 
and mapping of the ionosphere. Observations of the ionosphere 
from above by means of satellite-borne sounders offers the possi- 
bility of the rapid collection of synoptic ionospheric data on a global 
basis. Given data storage on the satellite and appropriate ground 
data transmission facilities, short-term ionospheric mapping may 
be possible. Results so far obtained from the Alouette satellite 
regarding the structure and variability of the topside ionosphere 
will be presented. 


Thermal conductivity and viscosity of simple fluids, J. V. 
Sengers, Intern. J. Heat Mass Transfer 8, 1103-1116 (1965). 

The behavior of the thermal conductivity and viscosity coefficient 
of simple fluids as a function of pressure, density and temperature 


is surveyed. Particular attention is given to the behavior of the 
transport coefficients in the critical region. In a previous paper [1] 
it was shown that the thermal conductivity coefficient exhibits a 
pronounced maximum in the critical region. Additional information 
confirming the existence of such a maximum is presented. 

Theoretical predictions of the transport coefficients in dense systems 
are often based on the theory of Enskog. A survey is presented 
showing to what extent the theory of Enskog describes the density 
dependence of the transport coefficients in monatomic fluids. 
Finally, the transport mechanism leading to the anomalous rise of 
the thermal conductivity in the critical region is indicated. 


International comparison on _ dielectric measurements, 
E. Rushton, G. Russell, B. W. Petley, H. E. Bussey, J. E. Gray, 
E. C. Bamberger, and D. Morris, JEEE: Trans. Instr. Meas. IM-13, 
No. 4, 305-311 (Dec. 1964). 

Three materials, namely fused silica, glass, and alumina, were 
selected for comparison based on known or expected homogeneity, 
isotropy, and stability. Measurements were made by the three 
international laboratories, both in the radio frequency range using 
capacitor type holders either with or without an air gap, and at 
microwave frequencies using either cavity resonance methods or 
transmission line impedance methods. The results show agree- 
ment among the three laboratories of better than one per cent for 


The large dimen- _ 


the real part of the permittivity. Also, the dispersion from 9 GHz 
to 1 MHz is reasonable. The agreement on loss tangent is less 
satisfactory and is being investigated. A second round of measure- 
ments is in progress and with the expectation that closer agreement 
will result. 


Lightning characteristics as derived from sferies, W. L. 
Taylor, (Proc. 3rd Intern. Conf. Atmospheric and Space Electricity, 
Montreus, Switzerland, May 1963), Book, Problems of Atmospheric 
and Space Electricity, ed. S. C. Coroniti, pp. 388-404 (Elsevier 
Publ. Corp., Amsterdam, The Netherlands, 1965). 

Lightning discharges emit electromagnetic signals, called atmos- 
pherics or sferics, that usually attain a spectral peak within the VLF 
region; i.e., 3 to 30 ke/s. Primary emphasis is placed on this band 
of frequencies for describing lightning discharge characteristics. 
The important problems and limitations in utilizing atmospherics 
are discussed with reference to particular research areas. These 
areas include: 1) measurements of the radiation field from light- 
ning discharges, 2) determination of VLF and ELF propagation 
characteristics, 3) observations of atmospherics for specialized 
purposes such as study of the thunderstorm mechanisms and detect- 
ing and tracking thunderstorms, and 4) radio interference from 
lightning discharges. 


New developments in dental materials: A world-wide survey, 
G. C. Paffenbarger, Intern. Dental J. 15, No. 3, 311-422 (Sept. 1965). 
General discussion of a symposium prepared by the moderator, 
G. C. Paffenbarger, chairman, Commission on’ Dental Materials, 
Instruments, Equipment and Therapeutics of the Federation 
Dentaire Internationale. 


MAGIC—A machine for automatic graphics interface to a 
computer, D. E. Rippy and D. C. Humphries, Proc. Fall Joint 
Computer Conf., Nov. 20-Dec. 2, 1965, Los Vegas, Nevada, pp. 
819-830 (Spartan Books, Inc., Wash., D.C. 1965). 

This report describes a machine which has been developed within 
the Computer Technology Section of the National Bureau of Stand- 
ards as a research tool for the investigation of man-machine com- 
munication techniques. This machine has been designated MAGIC 
(Machine for Automatic Graphics Interface to a Computer) and 
consists of large-diameter cathode-ray displays combined with a 
small, specially designed, programmable digital computer. It is 
intended to operate as a remote display station connected to a large 
ADP system via voice-quality communication lines. 

MAGIC is subdivided into two basic units: the processor unit and 
the display unit. The display unit consists of two CRT. displays 
termed primary and secondary. ‘The processing unit is also sub- 
divided into a control processor and four subordinate list processors. 
Three of these subordinate list processors operate directly on data 
displayed by the primary display, utilizing list processing techniques. 
Described is the software repertoire and hardware design of MAGIC, 
particularly that which pertains to the subordinate list processors. 
Conclusions are drawn with respect to the adequacy of the original 
system design. Future plans for MAGIC, and the proposed design 
and construction of a Model If MAGIC are discussed. 


The NBS time scale and its relation to other time scales, 
J. A. Barnes and R. C. Mockler, (Proc. XIV General Assembly of 
URSI, Tokyo, Japan, 1963), Book, Progress in Radio Science 1960- 
1963. Vol. 1. Radio Standards and Measurements, ed. R. W. 
Beatty, p. 41 (Elsevier Publ. Co., Amsterdam, The Netherlands, 
1965). 

Newman, Fey and Atkinson established an atomic time scale using 
the time pulses of W WV and based on the United States Frequency 
Standard (USFS) in Boulder, Colorado. The hyperfine structure 
separation was assumed to be 9192631770.00 ...cps. The rate of 
the WWV clock is determined by comparison with VLF trans- 
missions from Boulder which are in turn referred to the USFS. 
NBS atomic time was assigned to WWV time pulses beginning in 
October 1957. These assignments were compared with similar 
assignments by the U.S. Naval Observatory according to their 
A. 1 time scale. The NBS atomic time was made to agree with A. 1 
time on January 1, 1958. Since that date the two scales have 
diverged at an average rate of about 1 x 10~—!° sec/sec. 

More recently atomic clocks have been placed in operation at 
Boulder also based on the USFS. The precision of these clocks is 
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normally about 0.3 microseconds per day (3 x 10-'*). The ab- 
solute time differences between this clock at Boulder, the clock at 
WWVY, the clock at the U.S. Naval Observatory and the clock at 
the Loran-C master station at Cape Fear, North Carolina, were 
determined by carrying a portable clock between these various 
stations. The absolute time difference between these clocks was 
known to + 5y sec. at the time of the measurement, April 24, 1963. 
The results of this experiment yielded experimental values for the 
composite propagation and receiver delays between these stations. 
With these results the Loran-C time pulses are now used to con- 
tinuously compare the clock times of NBS-A and the U.S. Naval 
Observatory time A. 1. 


Annual report of the National Bureau of Standards, C. E. 
Moore, Astron. J. 70, No. 9, 637-639 (Nov. 1965). 

This report is furnished annually to the Astronomical Journal in 
order to present to astrophysicists the spectroscopic work that they 
are particularly interested in. It is published with the Observatory 
Reports and provides a brief summary of activity at the NBS. 


A review of zine oxide-eugenol type filling materials and 
cements, G. M. Brauer, Rev. Belge. Med. Dent.-Belg. Tijds. 
Tandheelk 20, No. 3, 323-364 (1965). 

\ review is given of the status of zinc oxide-eugenol type cements 
as well as modifications of these materials developed in this labora- 
tory and elsewhere. The mechanism of the setting reaction and the 
physical. histological and clinical properties of these cements are 
discussed. Zinc oxide-eugenol cement is decidedly superior to 
other materials in its palliative effect on the pulp and in patient 
comfort. Whereas. solubility is low and initial adaptation to the 
cavity walls-unusually good, its low strength and lack of resistance 
to abrasion confine its use to restorations not involving high stress 
or attrition. Properly formulated, it is the cement of choice for 
selected restorations. ' 
Recent investigations have explored the possibilities of improving 
the desirable physical properties. Incorporation of o-ethoxybenzoic 
acid has led to materials that appear promising as bases and cement- 
ing materials for inlays and crowns. Further studies should lead 
to stronger, more adhesive and more permanent restorative materials 
that are impermeable to moisture. 


Separation of pyrenediones by column chromatography, 
A. J. Fatiadi, J. Chromatog. 20, 319-324 (1965). 

The 1,6- and 1,8-pyrenedione were separated directly on a column 
of silica gel by means of glacial acetic acid as the eluant. The final 
purification of 1,6-pyrenedione was accomplished on a-‘column of 
activated alumina by use of benzene as the solvent. The IR and 
UV spectra of purified 1,6-, 1,8-, and 4.5-pyrenedione and of 1-oxo- 
6,7-phenalenedicarboxylic acid anhydride were recorded. 


The thermal degradation mechanism of polystyrene, L. A. 


Wall, S. Straus, J. H. Flynn, D. McIntyre, and R. Simha, J. Phys. 
Chem. 70, No. 1, 53-62 (Jan. 1966). 

New and extensive measurements of (1) the rates of volatilization 
and (2) the molecular weights as a function of the extent of pyrolysis 
are reported for various polystyrene samples, in particular, frac- 
tions spanning a broad range of molecular weights. The experi- 
mental observations are shown to be described quite well by the 
results of computer calculations based on the theory of free-radical 
degradation and a mechanism involving four elementary reactions — 
initiation, depropagation, intermolecular transfer, and termination 
by disproportionation. Certain characteristics of the volatilization 
rates are sensitive to intermolecular transfer, which is required 
together with pure end initiation in order to match theory with 
experiment. “Weak links” appear to be absent although intro- 
duction of a relatively slight amount of random scission, which would 
be important at large molecular weights, may improve the agree- 
ment. A tendency toward higher than theoretical weight-average 
degrees of polymerization during pyrolysis may be indicative of a 
small amount of termination by combination, diffusion-controlled 
rates, and a failure of the molecular weight relationship. 


Zone refiner with temperature control, 
Chem. 37, 1178-1180 (Aug. 1965). 

The sample tube is heated and cooled by fluids of controlled tem- 
peratures, circulating through successive chambers in a carrier 


F. Ordway, Anal. 


that slides along the rotating sample tube. Mixing and loss of the 
fluids are substantially eliminated by use of commercial shaft 
seals to separate the chambers. 


A simple method for making stereoscopic drawings, F. 
Ordway, Am. Mineralogist 50, 1496-1498 (1965). 

Stereoscopic pairs of drawings are made on a film-measuring de- 
vice that has a cursor movable along a graduated scale. One draw- 
ing is fixed to the scale, the other to the cursor. Points are plotted 
simultaneously on both by using carbon paper or a needle point, 
after setting the cursor for the required parallax distance. 


A system for accurate direct and alternating voltage meas- 
urements, F. L. Hermach, J. E. Griffin, and E. S. Williams, EEE 
Trans. Instr. Meas. IM-14, No. 4, 215-224 (Dec. 1965). 

A system for calibrating and using a group of de and ac standards 
has been developed at NBS to meet increasing needs for greater 
accuracy in the measurement of direct and rms alternating voltages 
at audio frequencies. With a group of saturated cells, a universal 
ratio set as a potentiometer, a simplified self-calibrating volt box, 
and a new differential-thermocouple ac-de comparator, an accuracy 
of 10 ppm can be obtainéd for direct voltage measurements and 40 
ppm for rms alternating voltage measurements between 20 and 
20,000 c/s, in terms of the direct volt and the ac-dc transfer standards 
maintained at NBS. 


Effect of surface roughness on emittance of non-metals, 
J. C. Richmond, J. Opt. Soc. Am. 56, No. 2, 253-254, (Feb. 1966). 
It has been observed experimentally, both at NBS and elsewhere, 
that the emittance of polished metals can be markedly increased 
by roughening the surface, by as much as a factor of 2 or 3. For 
non-metals and particularly white ceramic materials, on the other 
hand, the emittance appears to be essentially independent of sur- 
face roughness, at least for wavelengths below 7 or 8 microns. 
This apparent anomaly is explained on the basis of the differences 
in the optical properties of the two types of materials. 


Comparison of uncertainty ellipses calculated from two 
spectrophotometric colorimetry methods by an automatic- 
computer program, |. Mineroff, J. Opt. Soc. Am. 56, No. 2, 
230-237 (Feb. 1966). 

The parameters and equations for calculating uncertainty of chro- 
maticity coordinates in spectrophotometric tristimulus colorimetry 
have already been derived for the 10°-field standard-observer system. 
From the variability data for eleven optical filters obtained from two 
methods of spectrophotometric tristimulus colorimetry, standard- 
observer-system method and actual-observer method, two sets of 
chromaticity uncertainty ellipses are derived and intercompared. 
The similarity between corresponding ellipses of these two sets of 
ellipses leads to the conclusion that the system-derived uncertainty 
ellipses, based on the variances and covariances of the spectral 
tristimulus values of the large-field standard-observer system, can 
be used to estimate observer-derived chromaticity-uncertainty 
ellipses. : 
The equations which are the basis of the automatic-computer pro- 
gram and with which to compute uncertainty ellipses are given to 
aid anyone who wishes to arrange such a program in his laboratory. 
The fixed combinations of the spectral tristimulus values of their 
variances and covariances can be readily prepared by researchers 
who may wish to use desk calculators to compute uncertainty ellipses 
for actual-observer data or for standard-observer-system data, 
respectively. 


An international comparison of current-ratio standards at 
audio frequencies, B. L. Dunfee and W. J. M. Moore, JEEE 
Trans. Instr. Meas. IM-14, No. 4, 172-177 ( Dec. 1965). 

The results and analysis of an intercomparison between the National 
Research Council (NRC), Canada and the National Bureau of 
Standards (NBS), Washington, of two contrasting types of current 
ratio standards (current comparator and current transformer) are 
presented. The agreement achieved between the two laboratories 
when their respective designs and methods of measurement were 
quite different is emphasized. To this end, the basic theories under- 
lying the design and operation of the current comparator and current 
transformer are contrasted; the origin and significance of their 
respective errors are summarized: the (V+1 ) method used at NBS 
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to measure the errors of the multi-range transformer standards up 
to 10 kHz is contrasted with the N+ 1 method (one of several) 
used at NRC in calibrating the comparator standards. The com- 
parison circuit used at NRC to compare the respective standards 
of the two laboratories is described. The errors of the transformers 
as measured at NBS and NRC are presented and contrasted in 
both tabular and graphical forms. 


The design and performance of multirange current trans- 
former standards for audio frequencies, B. L. Dunfee, IEEE 
Trans. Instr. Meas. IM-14, No. 4, 190-204 ( Dec. 1965). 

The design, construction, and performance of two nearly identical, 
multi-range current transformers for operation 400 Hz to 16 kHz 
and with errors of only a few ppm are discussed. A general dis- 
cussion of the parameters that influence the performance of single- 
stage transformers and their consideration in effecting minimum 
errors throughout a wide frequency band is included. Construc- 
tion of the two standards having consecutive ratios 1/1 to 6/1 and 
operating at a rated secondary current of 5 amperes is described 
with the aid of line drawings. Circuit diagrams are presented and 
used to describe (1) the self-calibration circuit for measuring the 
errors of the 1/1 ratios and (2) the (N+ 1) circuit together with an 
auxiliary circuit for measuring the errors of ratios > 1. Balance 
equations are derived in an Appendix. Measured values are pres- 
ented in tabular form and include, in addition, the results obtained 
at the National Research Council Canada, during an international 
comparison and described in a companion paper. The effects of 
polarization at low frequencies is particularly emphasized by utiliz- 
ing the 400 Hz data together with information obtained in supple- 
mentary experiments; also, the role of winding and interwinding 
capacitances and the necessity for a more explicit definition of 
current ratio at the higher frequencies is high-lighted. 


High-frequency elastic moduli of simple fluids, R. Zwanzig 
and R. D. Mountain, J. Chem. Phys. 43, No. 12, 4464-4471 (Dec. 
15, 1965). 

This article presents a calculation of the infinite-frequency elastic 
moduli of monatomic fluids. When the intermolecular potential 
has the standard Lennard-Jones form, the elastic moduli are shown 
to be related to the pressure and internal energy of the fluid. Nu- 
merical values of the elastic moduli at various densities and tem- 
peratures are presented in tabular and graphical form. An identity 
is derived relating the shear and bulk moduli of any isotropic material 
in which particles interact by means of two-body central forces; 
this is a generalization of the familiar Cauchy identity occurring in 
the theory of elasticity of solids. The approach is based on analysis 
of the high-frequency limit of exact time-correlation function ex- 
pressions for shear and bulk viscosity. 


Photometry, L. E. Barbrow, Book, Encyclopedia of Physics, pp. 
516-519 (Reinhold Publ. Corp., New York, N.Y. 1966). 

The photometric quantities luminous intensity and flux, illuminance 
and luminance, and luminous reflectance, transmittance, directional 
reflectance factor, and directional transmittance factor are defined, 
the methods of measuring them are briefly described, and the factors 
which affect the accuracy of the measurements are discussed. 
Included is a table of the currently internationally agreed-upon 
symbols for the photometric and related radiometric quantities and 
units. 


Adhesive bonding of various materials to hard tooth tissues. 
V. The effect of a surface-active comonomer on adhesion 
to diverse substrates, R. L. Bowen, J. Dental Res. 44, No. 6, 
1369-1373 (Nov.-Dec. 1965). 

A surface-active comonomer or coupling agent (NPG-GMA, the 
addition-reaction product of N-phenylglycine and glycidyl meth- 
acrylate) has been found to improve significantly the bond strength 
between a direct-filling resin and fluorapatite. Since fluorapatite 
surfaces contain no organic material, this demonstrates that the 
collagen and other organic constituents of enamel and dentin are 
not essential to the improvement of bonding produced by the NPG- 
GMA. It does not follow, however, that the organic constituents 
play no role in the bonding. The objectives of the present in- 
vestigation were primarily to observe the effects of NPG-GMA on 
the bonding to unmineralized collagen and keratin and secondarily 


to learn more about NPG-GMA and its behavior in bringing about 
improved bonding between a methacrylate resin and hard tooth 
tissues. 


Chemical analysis, G. M. Brauer and G. M. Kline, Encyclopedia of 
Polymer and Technology 3, 633-665 (John Wiley & Sons Inc., New 
York, N.Y., 1965). 

The chemical analysis of polymeric materials is reviewed. Prepa- 
ration of the sample, physical tests, thermal degradation techniques 
and chemical procedures such as analysis for functional and end 
groups and characterization tests are discussed. Instrumental 
techniques that are employed in the analysis of polymers are de- 
scribed and systematic schemes for the analysis of high molecular 
weight compounds are outlined. 


Viscosity and density of boron trioxide, A. Napolitano, P. B. 
Macedo, and E. T. Hawkins, J. Am. Ceramic Soc. 48, No. 12, 613- 
616 (1965). 

Extensive and accurate viscosity values for molten BO; are reported 
over the continuous range from 20 poises (at 1400 °C) to 10" poises 
(at 318 °C). The measurements were made with a wide range 
rotating cylinder viscometer. The melts were pretreated by bub- 
bling super dry nitrogen through them to insure that the water con- 
tent was at a minimum. Densities were also determined between 
1400 and 400 °C. Above 800 °C the temperature dependence of 
the viscosity curve was Arrhenius, below this it had a smooth curva- 
ture showing no breaks. At the low temperatures, although its 
curvature decreased, the viscosity curve was not Arrhenius. 


Comparison of xenon and carbon ares as radiation sources 
for laboratory weathering of asphalts, K. G. Martin, P. G. 
Campbell, and J. R. Wright, Am. Soc. Testing Mater. Proc. 65, 
809-830 (1965). 

The xenon- and carbon-arc weatherometers were compared as 
radiation sources for the laboratory weathering of coating-grade 
asphalts using aluminum-backed asphalt films of 25 mils thickness 
and unsupported films of 25 microns (1 mil) thickness. The 25-mil 
films were exposed to the ASTM Daily Cycle A and degradation 
was measured by the spark-probe technique. The 25-micron films 
were irradiated at a controlled humidity level, but were not exposed 
to the cold water spray. Photooxidation was measured by infrared 
spectroscopy. 

Since the xenon source operated on a near continuous basis, it 
required less total time to cause failure of the 25-mil films than did 
the carbon-are machine, although the latter degraded asphalts about 
1.2 times faster per hour of exposure. The order of durability of 
asphalts, from a given source group, was reversed for the carbon 
arc compared to the xenon arc, but when asphalts of similar source 
were considered in groups to the overall correlation between the 
two weatherometer types was good. The spark-probe test method 
was found to influence greatly the durability results from both 
weatherometers. Thermal shock, alone, had little effect. Photo- 
oxidation rates of the 25-micron films were of the same order for 
xenon and solar sources, but a different order for carbon-are 
exposures. 


Proposed standard for measuring and reporting physical 
properties of optical materials, G. W. Cleek, F. J. P. Consitt, 
and W. D. Lawson, Infrared Phys. 5, 141-160 (Dec. 1965). 

For some time there has been evidence of a lack of data on the 
properties of optical materials, particularly as a function of tem- 


perature. This has been particularly apparent to designers of 
optical equipment which has to operate in unusual environments 
e.g., in space vehicles, or which has to withstand severe conditions 
during critical periods of operation e.g., during rocket launch. For 
instance refractive indices have in general been measured only 
around room temperature but optical instruments with complex 
lens systems may be required to function at either much higher or 
much lower temperatures. 

Another problem arises through different manufacturers of optical 
materials using different methods and specimens for measuring 
physical properties and presenting the data differently. This 
frequently makes it difficult for a user to compare directly optical 
materials from different sources and in general he does not wish 
to set up and carry out comparative measurements himself. Thus 
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the development of a joint set of standards on methods of measur- 
ing and reporting the properties of optical materials was considered 
highly desirable. 

With the needs of the optical designer in mind the authors have 
attempted to list specific properties of interest, outline standard 
methods of measuring these properties, and suggest standard 
methods of presenting the results. We have used for illustration, 
in Appendix C, the data on synthetic sapphire listed in the manu- 
facturer’s data sheets. 


Techniques for detection of high-altitude nuclear explo- 
sions, A. G. Jean, C. E. Hornback, R. H. Doherty, J. R. Winkelman, 
and R. T. Moore, Proc. IEEE 53, No. 12, 2088-2098 (Dec. 1965). 
The Central Radio Propagation Laboratory (CRPL) is conducting 
a study for the Advanced Research Projects Agency (ARPA) project 
VELA of ground based techniques for detecting geophysical effects 
produced by nuclear explosions at altitudes above about 20 km and 
in space. This activity is called HANDS for High Altitude Nuclear 
Detection Studies. A substantial amount of experimental data 
exists on effects produced by high-altitude explosions and theoretical 
predictions have been made of effects expected from explosions 
in space. The ability to record and identify these signatures is 
limited by knowledge of the natural background. The objective 
of the HANDS is to acquire data on natural background required 
in the design of future detection systems. The experimental 
facility for collecting background data at conventional and more 
rapid rates from a collection of colocated geophysical sensors is 
described. 


The analysis, synthesis, and description of biological images, 
L. E. Lipkin, W. C. Watt, and R. A. Kirsch, Ann. N.Y. Acad. Sci. 
128, No. 3, 984-1012 (Jan. 31, 1966). 

There may be distinguished two approaches to characterizing the 
properties of biological images: the statistical approach tradition- 
ally used for mechanized image-processing, and the articular ap- 
proach generally used among biologists. We discuss both 
approaches, but focus attention on the second, for we hold that an 
articular approach allows the expression of much that is effectively 
impossible to express in the form of numerical measurements; 
and we know of no inherent reason why a computer cannot deal 
with information of a nonquantitative nature. 

The articulation of biological images necessarily takes place on two 
fronts: the imposition of an articular structure on the image itself; 
and the expression of this construal in English sentences. We will 
study both kinds of articulation, and we will suggest how they may 
be brought together, in the form of linked pictorial and linguistic 
grammars. The computer system which this paper envisions will 
be able to analyze a presented image with respect to a pictorial 
grammar, and to formulate and accept descriptions of that image, 
in English sentences, with respect to a linguistic grammar. It 
will be able to present pictorial instances of English descriptions, 
and in other ways to respond to English directives. 


Multivibrator provides bidirectional output pulses, F. C. 
Ruegg, Electronics 38, No. 17, 87 (Aug. 8, 1965). 

A bipolar flip-flop with variable output levels has been designed. 
The circuit can be used to test the frequency response of many 
transistor circuits, and there are many applications in the field of 
logic circuitry. 


Precision of cubic lattice parameter measurements by the 
Kossel technique, H. Yakowitz, Trans. Symp. Electron Probe 
Analysis, Oct. 12-15, 1964, Wash., D.C., pp. 417-438 (John Wiley 
& Sons, Inc., New York, N.Y., 1966). 

The major factors affecting the precision of cubic lattice parameter 
measurement by means of the Kossel technique are considered. 
It is shown that three special cases can be solved rigorously for 
projection distortion effects and that obtaining these cases is not 
difficult. The fact that a lens ratio method should be employed to 
reduce measured lengths on the film to lattice parameter values is 
emphasized. General formulae for the lattice parameter and the 
geometric sensitivity are presented. Correction and minimization 
of errors resulting from refractive index effects and temperature 
variations is indicated to be possible. An example of the determina- 
tion of the lattice parameter of LiF is given. The indicated sensi- 


tivity of the method expresses the relative standard error of the mean | 
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lattice parameter of LiF, is given by =“ 
ati 


= 2.54 x 10-6, 
“ 


Status of linear relations among heats of transport, R. E. 
Howard and A. B. Lidiard, J. Chem. Phys. 43, No. 11, 4158-4160 
(Dec. 1, 1965). 

It is shown that previously derived linear relations among the heats 
of transport contain no physical information, but serve only to com-. 
plete the definition of these quantities in situations where the fluxes 
of the various components of the system are not linearly independent. 
Other defining equations are possible but physical predictions are 
unaffected by the choice so presented. This resolves the paradox 
of the apparent inconsistency of these linear relations with kinetic 
theories of heats of transport. 


Stress-strain properties of textile yarns subjected to rifle 
bullet impact, C. A. Fenstermaker and J. C. Smith, J. Appl. 
Polymer Sci. 1, 125-146 (1965). 

If a flexible filament is struck transversely by a rifle bullet, and the 
motion observed by flash photography. data can be obtained on 
the resulting strain distribution and on the time required to break 
the filament. The theory and experimental procedure for obtain- 
ing these data are described and test data are given for a sample 
of polyester yarn. These data include strain distributions at various 
times after impact and. at various impact velocities, and a plot of 
breaking time versus impact velocity. Data analysis shows that 
although creep and stress relaxation occurred within 50 w sec after 
impact, during the next 150 w sec the creep and relaxation effects 
were too small to be observed. At an impact velocity of 174 m/sec 
the strain at the point of impact changed from an initial value of 
1.08 percent to a 50 w sec value of 1.27 percent, and the correspond- 
ing stress from 15.26 g/tex to 14.51 g/tex. The yarn sample broke 
within 10 « sec when impacted at 620 m/sec velocity, but required 
100 w sec to break at an impact velocity of 420 m/sec. 


The effect of humidity on ionization measurements with 
cavity and free-air chambers, T. P. Loftus, B. Petree, and J. T. 
Weaver, J. Radiology 86, No. 1, 149 (Jan. 1966). 

Measurements are described that show the effect of humidity on 
the measurement of gamma radiation with cavity ionization cham- 
bers. It is found that the chamber current does not vary by more 
than 0.06 percent as the relative humidity is changed by more than 
50 percent. ‘A calculation, based on the assumption that the air 
and its moisture produce ionization independently, predicts a de- 
crease of about 0.4 percent for the measurement in moist air. The 
discrepancy is believed to be connected with the value used for the 
average energy required to produce an ion pair in moist air. 
Similar results were obtained with a free-air chamber exposed to 
60 kVcp x rays. A current increase of only 0.03 percent was ob- 
served when the relative humidity was changed from 48 to 13 per- 
cent, although the calculated increase is 0.3 percent. 


Other NBS Publications 


J. Res. NBS 70A3 (Phys and Chem.), No. 3 (May-June 1966), 
0 


Thermal relaxation and Brillouin scattering in liquids. R.D. Moun- 
tain (See above abstracts). 

Density of polyethylene crystals grown from solution. 
and E. Passaglia. 

Mechanical relaxation in polyethylene crystallized with various 
degrees of lamellar orientation. J.M.Crissman and E. Passaglia. 
(See above abstracts). 

Pentafluoropheny! alkyl and vinyl] ethers. 
Wall. 

Calibration of germanium resistors at low temperatures (2-20 
°*Kelvin). G. Cataland and H. H. Plumb. (See above abstracts). 

Dissociation pressure of aluminum carbide using a rotating Knudsen 
cell. E.R. Plante and C. H. Schreyer. 

Determination of intermolecular potential fungtions from macro- 
scopic measurements. M. Klein. 


G. M. Martir 


W. J. Pummer and L. A. 
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J. Res. NBS 70B2 (Math. and Math. Phys.), No. 2 (Apr.- 

June 1966), 75 cents. 

On a sequence of points of interest for numerical quadrature. S. 
Haber. (See above abstracts.) 

On abstract numerical integrations. J. J. Sopka. 

Treatment of outliers in samples of size three. F. J. Anscombe and 
B. A. Barron. (See above abstracts.) 

A note on contaminated samples of size three. T. A. Willke. 

Realizing the distance matrix of a graph. A. J. Goldman. 

Pairs of nonsingular matrices. K. Goldberg. 


Upper bounds for the determinant of a row stochastic matrix. K. 
Goldberg. 
Groups preserving ordering in vectors. K. Goldberg. 


A table of radiation characteristics for uniformly spaced optimum 
endfire arrays with equal sidelobes, M. T. Ma and D. C. Hyovalti, 
NBS Mono. 95, (Dec. 10, 1965), 45 cents. 

Quarterly radio noise data September, October, November 1964, 
W. Q. Crichlow, R. T. Disney, and M. A. Jenkins, NBS Tech. 
Note 18-24 (Feb. 10, 1966), 50 cents. (Formerly the Central 
Radio Propagation Laboratory of the National Bureau of Stand- 
ards, now Environmental Science Services Administration, 
Boulder, Colo.) 

Materials for PLACEBO V, W. C. Watt, NBS Tech. Note 281 (Jan. 
17, 1966), 50 cents. 

An atlas of solar flare effects in the ionosphere observed with a 
high-frequency Doppler technique, September 1960—December 
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